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VORTEX-LATTICE FORTRAN PROGRAM FOR ESTIMATING 
SUBSONIC AERODYNAMIC CHARACTERISTICS 
OF COMPLEX PLANFORMS 

By Richard J. Margason and John E. Lamar 
Langley Research Center 

SUMMARY 


A FORTRAN computer program has been developed for esUmating ^ 
aerodynamic characteristics of complex planforms. The program represen s Ittang 
Dlanforms with a vortex lattice. These complex planforms Include wings with variable- 
sweep outer panels, wings with several changes In dihedral angle across the span, wings 
with twist and/or camber , and a wing In conJuncUon with either a tail or ^ana . 
aerodynamic characteristics of interest are Uft and pitching moment for both the flat 
and/or twisted wing, drag-due-to-llft parameter, leading-edge thrust, leading-edge 
Uol Ltltuons M leaLg-edge thrust and suctton coefficients, distributions of several 
spa^ loading coefficients, distribution of lifting pressure coefficient damp.ng-tn-pitch 
parameter, damping-in-roU parameter, and lift coefficient due to pitch ra e. 

This paper is intended as a user’s guide for program appUcaUon and sample cases 
are Included to Illustrate most of the options available lor use in the 
included is a study of the effect of the vortex-latUce arrangement on some the com- 
puted aerodynamic characteristics along with some recommendations for specify g 
vortex-lattice arrangements for particular types of planforms. 

INTRODUCTION 

in recent years, some wings have become very complex because of the varied speed 
regimes in which they are required to operate. Such wings may have variable sweep, 
several changes in dihedral angle across the span, or even a variable dihedral angle near 
the wing tip. Computing procedures for predicting the aerodynamic characteristics of 
these Wings become very involved if an adequate representation of the planform is to b 
made tL problem becomes more Involved when the body or body and tail are include 
in the representation. In order to solve this problem for preliminary designs or for 
parametric evaluaUons, a computer program has been developed for esUmating 
dynamic characteristics of these complex planforms. 



In this FORTRAN computer program, the planform in steady subsonic flow is 
represented by a vortex lattice. Although this type of representation is not new (for 
example, refs. 1 to 12), the present program has several useful features that are not 
found together in other generally available programs of either the vortex -lattice or 
pressure-doublet type (refs. 13 to 15). 

The program uses a minimum of input data to describe relatively complex plan- 
forms. These planforms may be described by up to 24 line segments on a semispan. 
They may have an outboard variable-sweep panel or they may have several dihedral 
angles across the span. In addition, two planforms may be used together to represent a 
combination of wings and tails or wing, bodies, and tails. The analysis in the present 
paper has been extended to handle planforms in a sidewash field. These velocities occur 
when a planform has dihedral or when a second planform is placed at a different height 
from the first planform. 

The program described in the present paper was developed from a basic program 
written several years ago, which has had considerable use at the Langley Research 
Center. In recent years this basic program has also been used in industry. The results 
have shown good correlation with experimental data. 

The present paper is intended to serve both as a description of the program and as 
a user's guide for its application. This paper describes in detail the program input data 
(appendix A) and output data (appendix B) and provides examples and typical running times 
of various types of configurations which can be handled (appendix C) along with a 
FORTRAN program listing (appendix D). In addition, the results of parametric applica- 
tions of this program are presented to provide guidance in specifying vortex-lattice 
arrangements which can be expected to give acceptable results. 

SYMBOLS 


The geometric description of planforms is based on the body-axis system with the 
origin on the planform center line. (See fig. 1 for positive directions.) The planform is 
replaced by a vortex lattice which is in a wind-axis system with the origin in the planform 
plane of symmetry. (See sketch (d) in text for details.) The axis system by which the 
geometric influence of a given horseshoe vortex is computed is wind oriented and referred 
to the origin of that horseshoe vortex (fig. 1). The units used for the physical quantities 
defined in this paper are given both in the International System of Units (SI) and in the 
U.S. Customary Units. For the purpose of the computer program, the length dimension 
is arbitrary for a given case; angles associated with planform are always in degrees. 

The symbols used for input data in the computer program are described in appendix A. 

The symbols used in the description of the program are defined as follows: 
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A 

Bk 

b 

CD,i 


aspect ratio; listed as AR in computer program output 
element of boundary-condition matrix, 47 t% 


wing span, m (ft) 

, X Induced drag 
induced drag coefficient, ^ g^ — — 


Cd induced drag parameter based on Munk's far-field solution 

Cd induced drag parameter based on near-field solution 

Ct lift coefficient, L/q^Sref 


'L,T 




lift coefficient based on additional loading and actual planform area 

lift coefficient due to pitch rate, — t r, per rad 

g /qcre f 

\ 2U 

lift-curve slope, P®*" 


ro lli ng-moment coefficient. 


Rolling moment 
<looSref*^ 


Cm 


SCj 

damping-in-roll parameter, . . per rad 

i%) 


- . _ Pitching moment 

pitching-moment coefficient about Y-axis, — 

Qoo'^ref^ref 


aCm/^CL longitudinal stability parameter 


'mr 


damping-in-pitch parameter , 


ac 


m 


"l^ref 


2U 


, per rad 


Cn element of circulation term matrix, Tn/u 
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ACp incremental pressure coefficient, ^ower — ^upper _ ^ 

^ Qoo Qco 

Cc leading-edge suction coefficient, ^^^^ion 

QooSref 

Crj. leading-edge thrust coefficient, Leading-gdge thrust 

^oo^ref 

c chord, m (ft) 

c^v average chord, Sj/b, m (ft) 

C(. chord along left trailing leg of elemental panel, m (ft) 

Cd section induced drag coefficient based on near -field solution 

c^ section lift coefficient 

*^ref reference chord, m (ft) 

Cg section leading-edge suction coefficient 

c^ section leading-edge thrust coefficient 

dii section induced drag based on near-field solution, N/m (Ib/ft) 

F influence function which geometrically relates influence of single horseshoe 

vortex to a quantity which is proportional to velocity induced at a point, 
m-1 (ft-1) 

F sum of influence function F at a control point on wing caused by two sym- 

metrically located horseshoe vortices, one on left half of wing and one on 
right half of wing, m-1 (ft-1) 

Gn,k element of influence function matrix, F^^n^k - ^v,n,k 

L lift for entire wing, N (lb) 

I lift per unit length of span, l/(2s cos (p), N/m (Ib/ft) 

4 


I lift per unit length of vortex filament, N/m (Ib/ft) 

I lift generated along a finite length of vortex filament, N (lb) 

My pitching moment for entire wing about Y-axis, m-N (ft-lb) 

Moo free-stream Mach number 

my pitching moment about Y-axis due to lift developed on elemental panel, m-N 

(ft-lb) 

N maximum number of elemental panels on entire wing 

Nq number of elemental panels in a chordwise row 

Ng number of chordwise rows of elemental panels on wing semispan 

p roll rate, rad/sec; also, pressure, N/m2 (lb/ft2) 

q pitch rate about Y-axis, rad/sec 

q„o free-stream dynamic pressure, N/m2 (lb/ft2) 

Sref reference area, m2 (ft2) 

S,- actual planform area, m2 (ft2) 

s horseshoe semiwidth in plane of horseshoe vortex, m (ft) 

” ®ref^(^®n cos 

t section leading-edge thrust per unit span, N/m (Ib/ft) 

U free-stream velocity, m/sec (ft/sec) 

u backwash velocity, m/sec (ft/sec) 

V resultant velocity, m/sec (ft/sec) 

v sidewash velocity, m/sec (ft/sec) 
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w 

X,Y,Z 

X,Y,Z 

X,Y,Z 

x,y,z 

x,y 

x,y,z 

*c/4 

^3 c/4 
X' = x//3 

yep 

O! 

/3 

r 

y 

Ar 

V 


downwash velocity, m/sec (ft/sec) 
axis system of a given horseshoe vortex (see fig. 1) 
body -axis system for planform (see fig. 1) 
wind-axis system 

distance along X-, Y-, and Z-axis, respectively, m (ft) 
distance along X- and Y-axis, respectively, m (ft) 
distance along X-, Y-, and Z-axis, respectively, m (ft) 
midspan x-location of quarter-chord of elemental panel, m (ft) 
midspan x-location of three-quarter-chord of elemental panel, m (ft) 

fractional spanwise distance from root chord to center of pressure on left 
wing panel 

angle of attack, deg 

induced angle of attack, rad 

Prandtl-Glauert correction factor to account for effect of compressibility in 
subsonic flow, yl - 

vortex strength, m2/sec (ft2/sec) 

r 

nondimensional lift, — or 

bU 2b 

net vortex strength along left trailing leg of elemental panel, m2/sec 
(ft2/sec) 

nondimensional spanwise coordinate, y/(b/2) 
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p density, kg/m3 (slugs/ft3) 

<p dihedral angle, in Y-^ plane, deg 

A planform leading-edge sweep angle, in X-Y plane, deg 

quarter-chord sweep angle of elemental panel; because of the small angle 
assumption, also used as sweep angle of spanwise horseshoe vortex fila- 
ment, in X-Y plane, deg 

i//' = tan~l((tan \p')/0) 

Subscripts: 

a additional; or angle of attack 

B twist and/or camber at = 0 for chordwise row of elemental panels 

b twist and/or camber at Cj^ = 0 for elemental panel 

d desired 

i index for elemental panel in chordwise row 

j maximum number of elemental panels in chordwise row 

k index for control point 

I left half of wing 

lower lower surface 

n index for elemental panel on wing semispan 

o value taken at Cj^ = 0 

r right half of wing 

rad per radian angle of attack 
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s spanwise bound vortex element 

t chordwise bound vortex element 

tc twist and/or camber 

u backwash 

upper upper surface 

V sidewash 

w downwash 


BASIC CONCEPTS AND UMITATIONS 


The vortex-lattice method is used in this computer program to determine the aero- 
dynamic characteristics of planforms at subsonic speeds. This method is an extension 
of the finite step lifting-line method originally described in reference 16 and applied in 
reference 11. This method assumes steady, irrotational, inviscid, incompressible, 
attached flow. The effects of compressibility are represented by application of the 
Prandtl-Glauert similarity rule to modify the planform geometry. Potential flow theory 
in the form of the Biot-Savart law is used to represent disturbances created in the flow 
field by the lift distribution of the planform. It is assumed that in any plane parallel to 
the X-Z plane the vertical displacements which occur in the wing or wake are neglected, 
except when the boundary conditions at the control points are determined. 

The planform is divided into many elemental panels. Each panel is replaced by a 
horseshoe vortex. This horseshoe vortex has a vortex filament across the quarter-chord 
of the panel and two filaments streamwise, one on each side of the panel starting at the 
quarter-chord and trailing downstream in the free-stream direction to infinity. Figure 1 
shows a typical horseshoe-vortex representation of a planform. The boundary condition 
for each horseshoe vortex is satisfied by requiring the inclination of the fluid streamlines 
to match the angle of attack at the three-quarter-chord point of its elemental panel. The 
circulations required to satisfy this tangent flow boundary condition is then determined by 
solving a matrix equation. Then, the Kutta-JoukowskL theorem for lift from a vortex fila- 
ment is used to determine the lift from each elemental panel. These lift results are then 
summed appropriately to obtain lift, pitching moment, and other aerodynamic character- 
istics. A similar procedure called the near -field solution is used to compute leading- 
edge thrust, suction, and induced drag. This program ignores the effect of thickness. 
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The lifting-surface planform is represented for the computer program by a series 
of up to 24 straight segments which are positioned counterclockwise around the perimeter 
of the left half of the planform. Lateral symmetry is presumed. The lines start at the 
leading ec^e of the plane of symmetry, go along the leading edge to the left tip of the 
planform, return along the trailing edge, and end at the trailing edge of the plane of sym- 
metry. The preciseness of the x and y Cartesian coordinates and dihedral angles, 
given as input data, determines the accuracy of the planform representation. It is 
recommended that the planform coordinates listed in the second group of the geometry 
output data given in appendix B be plotted and examined after each computation to verify 
the accuracy of the planform representation. This check should be made before using the 
aerodynamic output data. 

There are a number of restrictions and limitations in the application of this com- 
puter program. These limitations are discussed in detail in the program description and 
are noted with the appropriate input variables in appendix A. For the convenience of the 
program user, a complete list of restrictions and limitations is presented. 

The restrictions in the first group apply to all planforms and are as follows: 

(1) A maximum of two planforms may be specified. For examples, see sample 
case 1 for one planform and sample case 2 for two planforms. 

(2) A maximum of 24 straight-line segments may be used to define the left half of a 
planform. The lateral separation of the ends of these lines can be critical when the 
horseshoe vortices are laid out by the computer program. For details of the lateral 
separation requirements, see pages 12 and 13. 

(3) The maximum number of horseshoe vortices on the left side of the configuration 
plane of symmetry is 120. When two planforms are specified, the sum total of the vor- 
tices in both is limited to 120. Within this limit, the number of horseshoe vortices in any 
chordwise row may vary from 1 to 20 and the number of chordwise rows may vary from 1 
to 50. For examples, see the sample cases in appendix C. 

The limitations that apply only to variable-sweep planforms are (1) there should 
always be a fixed-sweep panel between the root chord and the outboard variable-sweep 
panel, (2) the pivot cannot be canted from the vertical, and (3) no provisions have been 
made for handling dihedral in the geometry calculations for the variable-sweep panel or 
at the intersection of this panel with the fixed portion of the wing. 

The limitations that apply only to planforms which have nonzero dihedral angles or 
to two planforms which do not lie in the same plane are (1) the variation in local chord 
must be continuous from the tip chord to the root chord of each planform specified, (2) the 
number of horseshoe vortices in each chordwise row must be at least two, and (3) the 
number of horseshoe vortices must be constant over the semispan of each planform. 
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Restrictions on allowed values or codes for individual items of input data are described 
in appendix A. 

The calculations presented herein were made with a computer which used approxi- 
mately 15 decimal digits. For other computers with fewer significant digits, it may be 
necessary to use double precision for some of the calculations. In addition, it may be 
necessary to change some of the tolerances used in the program. These tolerances are 
mentioned in either the text or the program listing. 

PROGRAM DESCRIPTION 


This FORTRAN program is used to compute the following aerodynamic character- 
istics: Cl^ Cl at Q!=0, O' at Cl = 0, y^p, Cmo> ^Chi/^Cl, 

/ 2 ^ 

^D,ii/*^L j spanwise distribution of additional wing loading, spanwise distribution of wing 

loading due to twist and camber, and spanwise distribution of basic wing loading. In addi- 
tion, the following aerodynamic characteristics are computed for a specified lift coeffi- 
cient: the incremental pressure coefficient for each elemental panel, the spanwise dis- 
tribution of the combined basic and additional wing loadings, the configuration angle of 
attack, and the contribution of the major planform to lift coefficient and induced drag coef- 
ficient. At an angle of attack of 1 rad, the induced drag, leading-edge thrust, and suction 
coefficients are computed for the entire configuration by using a near-field solution. This 
program can also be used to compute or both Cl_ and Cjn (rotary derivatives). 

These quantities are described in detail in Part HI of the Program Description. 

The computation in this program for the aerodynamic characteristics is divided into 
three parts: Part I contains the required geometric calculations. Part n contains the cir- 
culation term calculations, and Part HI contains the final output terms, calculations, and 
answer listings. These three parts coincide with the three overlays in the FORTRAN 
computer program. The input data are described in detail in appendix A, and the output 
data are described in detail in appendix B. Several sample cases are given to illustrate 
the use of the program. Listings of the input data and computed results for these sample 
cases (appendix C), along with the FORTRAN computer program (appendix D) are given. 

PART I - GEOMETRY COMPUTATION 

The first part of the program is used to compute the geometric arrangement 
required to represent the planform by a system of horseshoe vortices and is divided into 
three sections. In Section 1, a description of the planform (group one of the input data in 
appendix A) is read into the computer. In Section 2, configuration details (group two of 
the input data) are read into the computer. In Section 3, the horseshoe vortex lattice is 


10 



laid out. When two planforms are used to describe a wing-body -tail configuration, each 
of these sections is repeated for the second planform. At the beginning of the geometry 
computation, a data card is read which describes the number of planforms (either 1 or 2), 
the number of configurations for which values are to be computed, and the reference val- 
ues for chord and area. 


Section 1. Reference Planform 

The planform is described by a series of straight lines which are projected onto the 
X-Y plane from the deflected planform as shown in figure 1 for a double-delta planform. 
The primary geometric data are the locations of the intersections of the perimeter lines, 
the dihedral angles, and an indication as to whether the lines are on a fixed or movable 
panel. The pivot location is also required for a variable-sweep planform. These data 
are described in group one of the input data (appendix A). For variable-sweep wings, the 
planform used for input should be the configuration with the movable panel in a position 
where the maximum number of lines required to form its perimeter are exposed. 

Section 2. Configuration Computations 

The particular configuration for which aerodynamic characteristics are sought is 
described by group two input data which are read here. These data include the following 
quantities: An appropriate configuration number, the number of horseshoe vortices 
chordwise, the nominal number of vortices spanwise, the Mach number, the particular lift 
coefficient at which the total span load distribution is desired, the sweep angle of the out- 
board panel for variable-sweep wings, a code to Indicate whether C^p should be com- 
puted, a code to indicate whether Cl^ and Cmq should be computed, and a code for 
each planform to indicate whether it is flat or whether it has twist and/or camber. The 
foregoing data are punched on one card for each configuration as described in appendix A. 

The number of horseshoe vortices used in each chordwise row (SCW) can be constant 
across the span or it can vary. If it is constant, simply indicate the number on the con- 
figuration card and this value will be used on each planform of the group one input. If it 
varies, use 0 and add the required input cards to define the table of values (TBLSCW (I)) 
described in appendix A. However, it is usually desirable to use a constant value the first 
time a planform is used in the program. For all but the most simple planforms, the pro- 
gram adds some extra rows of horseshoe vortices. (This is described in Part I, 

Section 3.) As a result, the number of chordwise rows actually laid out (SSW) is usually 
greater than the nominal number of rows (VIC) and it takes one run through the program 
to determine the exact number and location of the rows. 

The lift coefficient at which the total span load distribution (basic loading plus addi- 
tional loading) is desired will usually be between 0 and 1. However, if a value of 11 is 
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specified, an induced drag polar is computed. In this case, the program will provide val- 
ues of for 11 values of Cl from -0.1 to 1, as well as values of ACp and the 

total span load distribution at a Cl of 1. 

If a planform has twist and/or camber, additional data cards are required with the 
group two input data. These data are the local angles of attack in radians at the control 
points when the root-chord angle of attack is 0°. The control point of each elemental 
panel is at the midspan three-quarter-chord line. Generally, it is necessary to compute 
the vortex-lattice arrangement for the planform without twist and camber to determine 
the locations at which the local angles of attack are required. The order in which these 
data are provided is described in detail in appendix A. If a planform has no twist and/or 
camber, no additional cards are required for group two input twist data because the pro- 
gram will assign 0 for the values of the local angles of attack. If variations in the basic 
wing planform are desired for additional computer cases, they may be obtained by 
repeating only the group two input data with appropriate changes in any of the aforemen- 
tioned variables. 

For a variable-sweep planform, the angle which describes the sweep should be on 
the leading edge of the movable panel adjacent to the fixed portion. The intersection 
points and slopes for the planform in the desired position are then computed. For a fixed 
planform, the sweep-angle specification is not required because the program will use the 
unaltered basic planform. The planform breakpoints are checked to see whether any con- 
secutive pair in the spanwise direction is less than (b/2)/2000 apart. If this occurs, the 
points are adjusted to coincide with each other. The adjustment is necessary to avoid a 
poorly conditioned matrix which could result in biased results for the circulation terms. 
Although this adjustment is usually adequate for planforms with no dihedral, it may not be 
sufficient for wings having dihedral or for use of this program in computers which have 
fewer than 15 significant decimal digits. This problem is discussed in detail in Part I, 
Section 3. 

When two planforms are specified, the progi*am compares the spanwise location of 
the breakpoints on both planforms inboard of the tip of the planform with the shorter semi- 
span. If all the breakpoints coincide spanwise, no action is taken. However, if one plan- 
form has a breakpoint which does not occur on the other planform, an additional breakpoint 
is added to the other planform on its leading edge. This is done to force all trailing legs 
from the horseshoe vortices to occur at the same spanwise location, which keeps a trailing 
leg from one planform from passing close by a control point on the other planform and 
prevents unrealistic induced velocities at that control point. 

The program determines the planform area and span projected to the X-Y plane and 
uses these values to compute the average chord. Planforms which have a constant angle 
of dihedral from the root chord to the tip chord have an average chord which is independent 
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of dihedral angle. However, wings with more than one dihedral angle have an average 
chord which is dependent on the individual dihedral angles. 


Section 3. Horseshoe Vortex Lattice 

In this section, the procedure by which the horseshoe vortex lattice is laid out is 
described. The planform is divided chordwise and spanwise along the surface into trap- 
ezoi dally shaped elemental panels; one horseshoe vortex is assigned to represent each 
panel. The horseshoe vortices are similar to those described in references 11 and 16 
and are sketched in figure 2 for a typical panel. The horseshoe vortex is composed of 
three vortex lines: a bound vortex which is swept to coincide with the elemental-panel 
quarter-chord sweep angle in the plane of the wing and two trailing vortices which extend 
chordwise parallel to the free stream to infinity behind the wing. Figure 1 shows a typi- 
cal chordwise row of horseshoe vortices on an arbitrary planform. The nominal width of 
these horseshoe vortices is the total semispan in the plane of the wii^ divided by the var- 
iable VIC. (See appendix A.) 

The procedure for laying out the horseshoe vortices and the elemental panels is to 
begin at the left tip with the first chordwise row of vortices and then proceed toward the 
wing root. The actual spanwise locations of the chordwise rows of horseshoe vortices 
are adjusted so that there is always a trailing vortex filament at points where there are 
intersections of lines with breakpoints of the planform. This adjustment may cause the 
horseshoe vortex width to be narrower or wider than the nominal width. When a horse- 
shoe vortex has one trailing vortex filament which coincides with a breakpoint, the width 
of the horseshoe vortex may vary from 0.5 to 1.5 times the nominal width. When both 
trailing legs coincide with breakpoints, the width may vary from a maximum of 1.5 times 
the nominal width to a minimum width of (b/2)/2000, as described previously in Section 2. 
For wings with zero dihedral angles, good results can be expected for horseshoe vortices 
of these widths. However, for planforms having dihedral, the span loading results may be 
poor when narrow (less than 0.5 times the nominal width) horseshoe vortices exist. 

Hence, special care must be used in describing a planform with dihedral so that these 
narrow horseshoe vortices will not be used. The number of chordwise rows actually laid 
out is given by the variable SSW. 

In the chordwise direction, the horseshoe vortices are distributed uniformly and the 
number of vortices is given by either the variable SCW or TBLSCW (I). The maximum 
number of horseshoe vortices in the chordwise direction is 20 and in the spanwise direc- 
tion the maximum number is 50 on a semispan. However , the total number of horseshoe 
vortices (either the product of SCW and SSW or the sum of TBLSCW (I)) permitted by the 
program is 120 on a semispan. The exact number generated by the program depends on 
the values of VIC and SCW and on the details of the planform. As many as one additional 
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chordwise row of horseshoe vortices may be generated by the program at each breakpoint 
outboard of the root. Wings with dihedral must always have at least two horseshoe vor- 
tices chordwise; wings without dihedral may have only one. The most desirable spanwlse- 
to-chordwise horseshoe-vortex ratio is examined in that portion of the paper entitled 
"Effect of Vortex-Lattice Arrangement on Computed Aerodynamic Characteristics." 

The Prandtl-Glauert correction factor is applied to the x-coordinates and the tan- 
gents of the sweep angle of the horseshoe vortices at this point to account for compres- 
sibility effects. 

Parametric studies can be performed on optional features selected by repeating the 
group two input data. These parameters include Mach number, vortex-lattice arrange- 
ment, desired lift coefficient, distribution of twist and camber, and sweep angle for a 
variable-sweep planform. The optional features include the computation of the rotary 

derivatives C^ or Cj^ and Cm_. This computation is accomplished by repeating 
*'P Q, 4 

the information required by group two of the input data for each additional case. Any 
number of additional cases may be used for a given initial wing planform set. A few 
limitations for variable-sweep planforms which should be noted are (1) the pivot cannot 
be canted from the vertical, (2) no provisions have been made for handling dihedral in the 
geometry calculations for the variable-sweep panel or at the intersection of this panel 
with the fixed portion of the wing, and (3) there should always be a fixed-sweep panel 
between the root chord and the outboard variable-sweep panel. 

PART n - VORTEX-STRENGTH COMPUTATION 

The vortex lattice laid out in Part I is now used in place of the real wing to generate 
the same flow field as the wing and to determine the forces and moments acting on the 
real wing. To perform these functions, the flow must be constrained so that it does not 
pass through the vortex lattice at specified points. These points are called control points 
and are at the midspan three-quarter-chord line of each elemental panel. This flow con- 
straint is called the "no flow" condition and is equivalent to requiring that the flow be 
tangent to the real wing mean-camber surface. Simultaneous matching of the no flow 
condition at all the control points is used to compute the required vortex strengths. This 
can be conveniently expressed in matrix form as 

{c)=[g]-'(b} (1) 

where C^, Gji^k> are the elements of these matrices. 

The matrix represents the numerical values satisfying the boundary condi- 
tions which are presented in sketches (a) to (d) and equations (2) to (4). The traditional 
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representation for flat wings is shown in sketch (a) of a wing chord. 



w cos Q! - U sin a = 0 '' 

This boundary condition may be extended to represent wings with dihedral. This exten- 
sion is shown in sketch (b) , which is a view looking upstream toward the trailing edge of 
the left half of the wing span. 


V 



w cos a cos <Pi ~ ^ sin - U sin ot cos - 0 

A view looking upstream toward the trailing edge of the right half of the wing span 
(sketch (c)) presents a somewhat different combination of velocity vectors for the no flow 
condition from that just considered. 
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Sketch (c) 


w cos O' cos (pj. + V sin - U sin a cos = 0 (4 

In the geometry convention for this paper 

= <Pl = ~(Pr 

This relationship can be used to show that equations (2) and (3) are identical and have the 
form 


w cos a cos <p - V sin <^ - U sin a cos <p = 0 (E 

or, for small angles of attack, 

w - V tan cp ~ Uq! (g 

In the present formulation of a vortex lattice, the angle of attack in equation (5) refers to 
the flow at the control point for each elemental panel. The vortex lattice is located in a 
plane parallel to the free stream as shown in sketch (d). 
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Flow angle of attack at 
each control point 

Typical spanwise 
vortex filament 



Vortex-lattice traUing filaments 


Sketch (d) 


The 


downwash velocity for a particular horseshoe vortex can be expressed as 


w(x,y,z) = — F^(x^ 


where the downwash influence coefficient is 

(y tan - xQcos 4> 


^ (y ^ ^o,2,^(y2tan2^ + zZsecZ.//' - 2yx- tan V) - 2z cos 0 sin <p{y * x' tan i^') 

f(x- + s cos <!> tan i//')cos <(, tan i^' + (y -t- s cos »)cos -i- (z + s sin »)si^ 

[(x’ + SCOS (p tan i//’)^ + (y + s cos <pr + (z + s sin <i>) J 

(X- - s cos 0 tan 0')cos 0 tan 0' + (y - s cos 0)cos 0 + (z - s sin 0)sin ^ 

n 1/2 

.)2 + (z - s sin 


^{x' - s cos <}) tan + (y - s cos <^) 


y - s cos (p 


X’ - s cos (p tan 4/^ 


(y - s cos (p)^ + (z - s sin p)^ | - s cos 0 tan + (y - s cos + (z - s sin p) 

X' + s cos p tan p' 


y + s cos 


(y + s cos p)^ + (z + s sin p)^ ^ g tan p 

and the sidewash velocity can be expressed as 


')2 + (y + s cos 0)^ + (z + s sin 0)^] 


where the sidewash influence coefficient is 




X* sin 0-2 cos <p tan 


jlj r 7- 

(X') + (y sin ^>)2 + cos20(y2tan2,^- + z2sec2^ - 2yx' tan ^j,') - 2z cos ^ sin 0(y + x' tan ^') 


j (x> + s c os tan i//^)cos <A tan (y + s cos 0)cos 0 + (z -<• s sin 0)sin d> 

[(x' + s cos 0 tan + (y + s cos 0)2 + (z + s sin 0)2]^^^ 

_ ( xj - s cos 0 tan 0')cos 0 tan 0' + (y - s cos 0)cos 0 + (z - s sin 0)sin 0^ 

[(x’ - s cos 0 tan 0’)^ + (y - s cos 0)2 + (z - s sin 0)2]^^^ 

+ ^ — ^ sin. 0 J ^ X' - s cos 0 tan 0’ | 

(y - s cos 0)2 + (z - s sin 0)2 S (7^ 7 ~ :;TT79/ 

[(x - s cos 0 tan 00^ + (y - s cos 0)2 + (z - s sin 0)2] ^ 


z + s sin 


x^ + s cos 0 tan 0^ 


(y + s cos 0)2 + (z + s sin 0)2 S I ~ 

[(x + s cos 0 tan 0 ')^ + (y + s cos 0)2 + (z + s sin 0)2]^^^ 

Then, by using equations (7) and (9) equation (6) can be rewritten as 
- Fy tan = Ua 

For a vortex lattice of N elements, equation (11) can be expressed for a particulai 
trol point by 


n=l 


^w,n " ^v,n 


For symmetrical aerodynamic loading on each half of the wing, equation (12) may be 
expressed as 

N/2 

7 (^w,n - Fy,n tan = Ana 

n=l 


where 


and 


w,n 


= F. 


w,n(x',y,z,s,i//',0)jgj^ + ^w,N+l-n^*'»y’^>s>>^’»<^)right 

panel panel 


+ F 


panel 


( 10 ) 


( 11 ) 

con- 

(12) 


(13) 


(14) 
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^v,n = Fy jj(x',y,z,s,;//',0)j^j^ 

panel 


(15) 


Figure 1 shows the locations of elemental panels n and (N + 1 - n). The matrix which 
is solved by the program is then 




^w,n,k " ^v,n,k tan 



(16) 


where Ok describes the local angle of attack in radians at the control point. For the 
first solution, is that angle of attack due to twist and camber when the root-chord 
angle of attack is zero; for the second solution, the angle of attack is 1 rad for all 
the control points. 

As previously mentioned, this program can be used to compute the rotary stability 
derivatives C^p, Cniq, and • This computation is accompUshed by following the 

method outlined in reference 17 where the values of the boundary conditions of the second 
solution are changed to an equivalent quasi-steady-state rolling or pitching motion. For 
steady-state rolling at zero angle of attack, the boundary conditions lead to a linear twist 
whose angle variation across the span is 



(17) 


For this computation, if the tip angle pb/2U is specified to be 5°, then equation (17) can 
be written as 


^ (2) = (18) 

' 2V\^/2j 180\b/2y 

For pitching motion, the Y-axis is the center of rotation. It is recommended that the 
perimeter points be specified so that the Y-axis coincides with either the center of gravity 
or the wing quarter -chord. For steady pitching motion, the boundary conditions lead to a 
parabolic camber as can be seen from 


«k(2) = 


-qx 

u 


-8Z 

9x 


(19) 


Specifying that 

q. ^ ^ (20) 

U 180 


leads to 


0^(2) = 


-Sttx 

180 


( 21 ) 
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If any of the rotary derivatives are to be computed, the program assigns zero values for 
the 0 ^( 1 ) terms and the appropriate boundary condition values for the q![j( 2) terms. 

In addition to solving for the circulation, solutions for section induced drag and 
leading-edge thrust are made at this point in the program by using a near-field approach. 
A detailed description of this implementation is given in Part HI, Section 3. 

PART ni - AERODYNAMIC COMPUTATION 

The circulation terms Pn/U computed in Part n are used in this part of the pro- 
gram to compute the lift and pitching-moment data for planforms with dihedral. A sim- 
plified procedure is used for zero-dihedral planforms. Then, the final form of the out- 
put data is obtained and printed for both planforms. 

The procedure described in Section 1 is used for planforms with dihedral and for 
wing-tail planforms where the planforms are not at the same elevation. A special treat- 
ment is needed for both types of planforms because there are local sidewash and back- 
wash velocities in addition to the free-stream velocity. The interaction of these velocity 
components with the spanwise bound vortex provides an additional lift force and the inter- 
action of the sidewash with the chordwise bound vortex (that portion of the horseshoe vor- 
tex trailing leg ahead of the wing trailing edge) results in another and new lift force. 
Because of the computation procedure used in Section 1, these types of planforms must 
have a continuous variation in local chord from the wing tip to the wing root. As a result, 
streamwise perimeter edges can only be used at the wing tip or tip of the tail for these 
planforms. 


Section 1. Lift and Moment Using Entire Horseshoe Vortex 

The lift, pitching-moment, and rolling-moment output data for planforms which have 
a nonzero dihedral angle over any portion of the planform or for two planforms at differ- 
ent elevations are computed here by using the local sidewash and backwash velocities in 
addition to the free-stream velocity. 

The procedure described herein for computing lift and pitching-moment data is per- 
formed twice: first, for the circulation terms due to twist and camber and, second, for 
the circulation terms due to an angle of attack of 1 rad. The lift, pitching-moment, and 
spanwise center-of-pressure data are computed for all elemental panels in a particular 
chordwise row; the procedure is then repeated for each chordwise row until the entire 
left half of the wing has been taken into account. For each elemental panel, the lift devel- 
oped along the left chordwise bound vortex is computed first and then the lift along the 
spanwise bound vortex is computed. The Kutta-Joukowski theorem for lift per unit length 
of a vortex filament is used to compute lift for wings with dihedral and is given by the 
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following equation: 


I = pVr (22) 

The circulation and velocity values used in equation (22) by this computer program are 
described in the discussion that follows. 

The lift developed along the chordwise bound vortices in a chordwise row of horse- 
■shoe vortices varies from leading edge to trailing edge of the wing because of the longi- 
tudinal variation of both the sidewash velocity and the local value of vortex strength. In 
figure 3, it can be seen that there is no circulation along the chordwise bound vortex from 
the leading edge of the wing to the quarter-chord of the first elemental panel. As a 
result, no lift can be generated here. On the chordwise bound vortex from the quarter- 
chord of the first elemental panel to the quarter-chord of the second elemental panel, 
there is a constant value of circulation and a varying value of sidewash velocity. A spe- 
cial case occurs for the first elemental panel at the left wing tip; there the value of cir- 
culation just equals that of the first elemental panel of the first chordwise row of horse- 
shoe vortices. Inboard from the tip, this chordwise bound vortex lies between two chord- 
wise rows of horseshoe vortices, and its circulation is equal to the difference between the 
circulations of the first elemental panel of each row. The sidewash velocity used is the 
one computed at the three -quarter -chord on the left chordwise bound vortex of the first 
elemental panel. 

The next lift to be computed is that developed along the chordwise bound vortex 
between the quarter-chord of the second elemental panel and the quarter-chord of the 
third elemental panel. This lift is computed in a manner similar to that of the first 
horseshoe vortex but there are differences and these are now explained. At the left wing 
tip, the sum of the circulation values of the first two elemental panels is used. Inboard 
from the tip between two chordwise rows of horseshoe vortices, the circulation is equal 
to the sum of the difference between the circulations of the first elemental panel of each 
row and the difference between the circulations of the second elemental panel of each row. 
The sidewash velocity used is the one computed at the three-quarter-chord on the left 
chordwise bound vortex of the second elemental panel. 

This procedure continues through the last elemental panel in a chordwise row. 
However, the final chordwise bound vortex extends from the quarter-chord of the last 
elemental panel to the trailing edge of the wing so that its length is equal to only three- 
quarters of the length of the other chordwise bound vortices in the same chordwise row of 
horseshoe vortices. The sidewash velocity described in the foregoing procedure is given 
by the following equation: 
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( 23 ) 


N/2 

X. - _L Y p 

U 47t 4 U 
n=l 


Horseshoe vortex filaments or their extensions which go through the point at which 
the velocity is being computed are eliminated in the computer program from equation (23) 
since a line vortex filament cannot induce a velocity on itself. The lift generated along 
an elemental length of chordwise bound vortex divided by free-stream dynamic pressure 
and reference wing area is given by 

_ A _ = _2_AIcc^ (24) 

qSref S^ef U U 


where AT is the local value of circulation as described in the preceding paragraph and 
Cc is the chord or elemental length of the chordwise bound vortex. No lift is computed 
along the chordwise bound vortex at the root because the sidewash velocity is zero for 
symmetric loading and geometry. 

The lift along the spanwise bound vortex depends on the values of free-stream, 
backwash, and sidewash velocities and on the circulation at the elemental panel. The 
sidewash velocity is given by equation (23) and the backwash velocity is computed from 


U 


N/2 
J_ Y £n 
Av ^ V 


Fu,n 


n=l 


(25) 


where 


Fu,n = Fu,n(x',y,z,s,;^',c^)jgjt n- 

panel panel 

and the backwash influence coefficient is 




z COS <t> - y sin 


(x*)^ + (y sin (p)^ + cos2<^(y2tan2i^ + z^sec^\p - 2yx’ tan i//’) - 2z cos <f> sin (f>(y + tan ;//’) 


(x' + s cos (p tan i^')cos <p tan (y + s cos <^)cos (p + {z + s sin 0)sin <p 


[(x' + s cos <^> tan + (y + s cos <p)^ + (z + s sin 0)^] 


1/2 


(x' - s cos 0 tan 0')cos 0 tan 0' + (y - s cos 0)cos 0 » (z - s sin 0)sln 0 1 
[(x' - s cos 0 tan 0’)^ + (y - s cos 0)^ + (z - s sin 0)^] ^ 


(27) 
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Equations (8), (10), and (27) represent an extension of the original formulation by Glauert 
(ref. 16) for rectangular horseshoe vortices, the later formulation by Campbell (ref. 11) 
for a spanwise vorticity filament with sweep, and the recent formulation by Blackwell 
(ref. 12) for a rectangular horseshoe vortex with dihedral. In contrast, the present equa- 
tions represent a subset of the formulation by Rubbert (ref. 3) in that the trailing legs are 
constrained to the free-stream direction. 

A spanwise bound vortex filament is shown in figure 4 and the lift generated along 
this vortex filament comes from both the total axial velocity interacting with the compo- 
nent of the vortex filament parallel to the Y-axis (2s cos <p) and the sidewash interacting 
with the component of the vortex filament parallel to the X-axis (2s tan \J/ cos <p). The 
expression for this lift divided by free-stream dynamic pressure and reference area is 


^s 

<looSref 


2 r 


^ref 


U 


(2s) 


1 - — I + — tan 
U U ^ 


cos 0 


(28) 


The contribution of the lift of the elemental panel to pitching moment is given by 
<loo^refCref ^loo^ref ^ref q^oSyef c^ef 


(29) 


To get the total wing lift and pitching-moment coefficients, these terms are summed 
over all the elemental panels which represent the wing in the following manner: 


(30) 


My 

‘loo^ref^^ref 


= 2 



(31) 


There are two values for each of these quantities; one for the surface loading due to twist 
and camber and the other for the surface loading at 1 rad angle of attack. From these 
quantities, four output terms are obtained. The lift-curve slope per radian is the value 
given by equation (30) (i.e., the lift coefficient at 1 rad angle of attack). The Uft-curve 
slope per degree is 


C 


L« = 



(32) 
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The 


longitudinal stabiUty parameter about the origin of the X-axis for the wing is given by 


Mt 


_ 

8Cl 


Woo^ref^refy 


(33) 




The pitching moment at zero lift is 




9Cni/ L 


(34) 


'm 


° vloo^ref^^ref/ 


The center of pressure in a 
expression: 


'tc 'tc 

spanwise direction is computed from the following 


N/2 

I 

n=l 




(35) 


cp 


1/ L 


2\jloo®ref/^\2 


The span-load coefficients are obtained from the Utt along the spanwise and chord- 
wise bound vortices of each horseshoe vortex. Before converting the lift expressions to 
span-load coefficients, a few basic definitions should be emphasized. The l‘ft “ <=<1“^*- 
Uons (24) and (28) is lift In units of force developed over a span equal to the width of a 
horseshoe vortex. Therefore, lift per unit length of span is 


I = 


I 


(36) 


2s cos 4> 

The span-load coefficient for an elemental panel is developed as follows: 

I 


CiC 


qooC 


i 


®ref 


(37) 


CL^av <^L^av \qooSref/CL2Sn cos <t>c^y 


where 

St 

cav=-g- 


( 38 ) 
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and 


T = 


^ref 

cos 


so that 


f T 

^L*^av ^oo^ref ^L 


(39) 


(40) 


At a particular spanwise location, each of these lifts are summed chordwise and converted 
to span-load coefficients by the following equations: For lift along the spanwise bound 
vortex filament, 



(41) 


For lift along the chordwise bound vortex filament, 



(42) 


Figure 5 shows the spanwise distribution of the span-load coefficients obtained from equa- 
tions (41) and (42) for a wing with dihedral. The results of these equations must now be 
combined to get the final distribution. It is assumed that the span-load coefficient should 
be zero at the wing tip, a result which cannot be obtained by direct combination of the 
results of equations (41) and (42). Since the vortex-lattice procedure is a finite approxi- 
mation for the continuous variation of circulation across the wing span, each value of cir- 
culation represents the average value over the width of one horseshoe vortex. For this 
calculation, it is assumed that the circulation terms or span-load terms are correct only 
at the center of each row of horseshoe vortices. The lift along the spanwise bound vorti- 
ces is computed here and is used directly; whereas, the lift along the chordwise bound 
vortices is interpolated linearly to determine its value at the midpoint of each row. 

These two values of lift are then combined as illustrated in figure 5 to give the final span- 
wise distribution of span-load coefficients. 

In order to determine the damping-in-roll parameter of wings with dihedral, the lift 
distribution which results from the antisymmetrical span loading must be combined with 
the appropriate spanwise moment arm. This combination can be expressed as 
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c, = 


^ flcoSrefb 


N/2 N/2 


n=l 


n=l 


(43) 


and, thus, 


c 

/pb\ 57T/180 

^\2uj 


(44) 


Section 2. Lift and Pitching and Rolling Moments Using Only Spanwise Filament 

of Horseshoe Vortex 

The computation of the lift, pitching-moment, and rolling-moment output data for 
wings which have no dihedral over any portion of the wing is described in this section. 
All the lift is generated by the free-stream velocity crossii^ the spanwise vortex fila- 
ment since there will be no sidewash or backwash velocities. For a single elemental 
panel, the lift per unit length of vorticity is 


T = pUr cos 

Since the length of vorticity is 2s/cos \p, the resultant lift is given by 

i = T 

cos \}/ 

Then, the lift per unit of span is defined by 




and is nondimensionalized in the following form for later use as 


I _ 2 r 

^oo^^av ^av U 
For a chordwise row 


j 



(45) 


(46) 


(47) 


(48) 


( 49 ) 
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The total lift coefficient is obtained by integrating the lift over the span as given by 



or approximately by 
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L " 


8 

Sref 


N/2 


A U 




n=l 


( 50 ) 


(51) 


The lift-curve slope per radian is obtained from a lift coefficient based on the circulation 
terms obtained at 1 rad angle of attack. 

The longitudinal stability about Y-axis is given by 


8Cm 


1 

^ref 


N/2 

Y ra,n. 

2 -jrxs,nsn 
n=l 


n/2 



n=l 


The pitching moment at zero Lift is 


(52) 


N/2 

Cm 8 y ^tc,n 

“°"CrefSref ^ U 

n=l ^ 

The center of pressure in a spanwise direction is 


(53) 
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yrn ~ 


1 n=l 


y^y s 
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V 

Z -t^"^ 
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(54) 
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The span-load coefficient is 


^L^av 


j 



N/2 


r. 


u 



n=l 


( 55 ) 


The same procedure used to compute the damping-in-roll parameter for wings with 
dihedral can be used to compute C^p for zero-dihedral wing planforms except that the 
contribution of the chordwise bound vortex is eliminated. Thus, equation (43) becomes 




‘loo^ref^ 


iN/2 

^ 2(l)js.n2sn 


n=l 


and likewise 


(56) 


'h 57T/180 


(57) 


Section 3. Output Data Preparation 

This section of the program is used to compute the last portion of the data listed in 
the final output. These data include the damping-in-pitch parameter, the lift coefficient 
due to pitch rate, the induced drag parameter, the angle of attack for zero lift, the angle 
of attack for the desired lift coefficient, the basic span load distribution, and the addi- 
tional span load distribution. 

The pitch derivatives can be computed by using the vortex strengths obtained with 
the boundary condition values which represent a constant pitching motion. These vortex 
strengths are employed to compute a^iid Cm which, in turn, are used as follows: 


and 


C = ^ ^m 

/qc\ " 57 t Cref 
180 2 


C 

Lq 57T ^ref 

'y2uj 180 2 


(58) 


( 59 ) 
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In this paper, induced drag parameters are computed by both far-field and near- 
field methods. The far-field method is based on the lifting-line concepts employed in 
the Treffetz plane by Munk and the induced drag parameter thereby obtained can be 
expressed mathematically as 






( 60 ) 


This equation has been reformulated by Multhopp using, in part, his quadrature formula 
and is programed here in the form presented by equation (146) in reference 18. Equa- 
tion (60) can give good results for wings without dihedral but should be used only as a 
guide for wings with dihedral, since no vertical displacement of the span loadings is 
taken into account. For wings having dihedral, a method such as that developed in ref- 
erence 19 or the near-field method should be used to compute the induced drag. Even for 
wings without dihedral, good results can only be expected for the far-field method when a 
large number of chordwise rows of horseshoe vortices are specified since the interpo- 
lating procedure chosen to represent the variation of y with sin-lrj was a linear 
curve fit between consecutive pairs of data points. This curve fit requires that a suffi- 
cient number of data points be available near the wing tip where the gradient of the 
y - sin'^Tj curve is the greatest. 

The near -field computation for the induced drag is based on combining for each 
elemental panel the lift and leading-edge thrust as follows: 


L (61) 

^loo ^loo ^loo 

where the lift per unit of span Z/q^ is computed by equation (48) for planforms without 
dihedral and by equations (24) and (28) for planforms with dihedral. The leading-edge 
thrust per unit of span is computed by using the Kutta-Joukov^id theorem where the 
induced and free-stream velocity components parallel to the Y-Z plane interact with the 
spanwise bound vortex filament as follows; 




- — tan (b - a 
U 



(62) 


There is no contribution of the chordwise bound vortex filaments to the leading-edge 
thrust. In contrast, however, there is a contribution of the lift due to the chordwise 
bound vortex filament included in the induced drag term. (See eqs. (6) and (24).) It 
should be noted that this equation is evaluated at an angle of attack of 1 rad and that the 
circulation used is the one due to the additional loading only. 
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These results are then summed along each chordAvise row to get the following sec- 
tion leading-edge thrust: 

2b 2b 4Uoo 


i=ix^°°/i 


From equation (63) the section suction coefficient is computed as 


CgC /C^c\ 


I cos A 


2b V2b,^ 

Then, the section induced drag for a chordwise row of horseshoe vortices is 
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Finally, the near -field solution for the induced drag parameter is 
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(66) 


In addition, the leading-edge thrust and suction coefficients are computed similarly as 
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= — / it 2sk cos 0k 
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and 
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( 68 ) 


The angle of attack for zero lift is computed by 


Ob = - 


^L,tc 

^L, 


0( 


(69) 


The angle of attack required for the additional loading and basic loading combined to pro- 
duce the input value of the desired lift coefficient is 


30 


(70) 


“d = 



+ «b 


The basic load due to twist and/or camber is the load on the wing when the lift 
coefficient is zero. This load is obtained from the values of c^c^Cav for each elemen- 
tal panel as follows; 



Equation (71) is then summed for each chordwise row for the span load distribution of 
basic load to give 



The span load distribution at the input value of desired lift coefficient is 



In addition, the span load distribution and local lift-coefficient ratio 

c^^Cl t listed where the lift coefficients are based on the lift due only to additional 
loading and the total lift coefficient Cl is based on the true planform area Sj. Also 
listed is the distribution of local chord ratio c/ c^v 

The incremental pressure coefficient is defined as 
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p,n 
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% 


Since the pressure is assumed to be uniform over an elemental panel, 
(^/c)n 


AC 


P,n--^ 


(74) 


(75) 


which is used in the program. For planforms without dihedral, equation (75) can be 
expressed as 
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EFFECT OF VORTEX- LATTICE ARRANGEMENT ON COMPUTED 
AERODYNAMIC CHARACTERISTICS 


Several sets of lifting-surface planforms have been investigated to determine the 
effect of the vortex-lattice arrangement on the computed aerodynamic characteristics. 
The first four sets of planforms had two prescribed leading-edge sweep angles in com- 
bination with three different taper ratios for aspect ratios of 2, 4.5, and 7. Calculated 
results for these planforms show that for different vortex-lattice arrangements, smaller 
variations of y^p and are produced than of and 

Cd u/Cl^- The variation of y^p with vortex-lattice arrangement is presented for 
unswept wings of taper ratio 1.0 in figure 6. These data indicate that increasing Ng 
leads toward converging results for y^^p for all N^. 

The variations of Cj__^ and Cj) with vortex- 

lattice arrangement are presented in figure 7 for unswept planforms with a taper ratio 
of 1.0 and in figures 8 to 10 for planforms with a leading-edge sweep angle of 45° and 
taper ratios of 1.0, 0.5, and 0, respectively. These data indicate the following conclu- 
sions. A spanwise increase in the number of chordwise rows of horseshoe vortices Ng 
leads to converging answers. For these simple planforms, the Ng required for con- 
vergence of to a particular value is sufficient for convergence of ^Chi/^Cl and 

^D,i/^L^ and should be 20 or larger. Also, the computed values of Cl^, 9Cm/®CL, 
and Cj) in most instances have a definite dependence upon Ng. In particular, 

Nc controls the asymptotic levels that these aerodynamic characteristics attain with 
varying Ns- These asymptotic levels approach a converged result when Ng is 
increased. Differences between asymptotic levels which occur for consecutive N^ val- 
ues decrease with increasing Ng and the largest difference in asymptotic levels is 
obtained by increasing Nc from 1 to 2. Therefore, an Nq value of 2 should be the 
minimum used. Higher values of Nq have little effect on Cj^ ; however, increasing 
Nc to 4 or more can provide additional improvement in ^Chi/^Cl and la 

contrast, the calculated results indicate that Nc has little effect on The 

asymptotic levels of Cj) and Cp) when Ns is greater than 20 can be 

compared with those of 1 /ttA. This comparison shows that converges to a 

value greater than I/ttA, as expected, whereas converges in a less uniform 

manner to a value less than 1/ nA. 
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Since is computed by using equations (65) and (66) which are based on 

ct and c;, these results indicate that ct may be overpredicted. However, a compari- 
son can be made in figure 11 between the distribution of section thrust computed for an 
A = 4 delta wing by the vortex-lattice and Wagner's (ref. 14) methods. It can be seen 
that the resulting magnitudes predicted by the two different methods compare closely in 
general shape and lead to comparable overall thrust results. From additional computer 
studies it has been found that the Nc = 10 and Ng = 12 pattern used for the results 
shown in figure 11 also provides reasonable results for other delta wings. The large 
number of chordwise stations is necessary on such wings so that the effect of the induced 
camber loading can be properly taken into account. Although the correct thrust coeffi- 
cient can be obtained from the far -field induced drag^and lift-curve slope directly, only 

by finding the appropriate combination of Nc and Ns will the induced-drag results be 

the same for both methods. This check provides a method by which the correct distribu- 
tion of section thrust can be obtained. The results presented in figures 7 to 10 show how 
difficult it is to make this check even for some simple planforms. 

To determine the effect of vortex-lattice arrangement on C^p, Cmq» and Cl^, 
additional computer studies were made with a cropped double-delta planform having an 
inboard leading-edge sweep angle of 83°, an outboard leading-edge sweep angle of 62°, 
and an aspect ratio of 1.49. Results of these studies showed two trends. :^r estimating 
Cl , a large value of Ng is desired with at least two horseshoe vortices (Nc) in each 
row. For estimating Cmq and CLq, a large value of Nc (8 or more) is desirable 

with a nominal value of Ng of 8 or 10. 

A final set of computer studies were made with the wing-body-tail configuration 
illustrated in sample cases 2, 3, and 4. The aerodynamic characteristics were computed 
ior this complex configuration by using 22 different vortex-lattice arrangements which 
had a total number of vortices on a semispan ranging from 17 to 120. Results showed 

very Uttle variation of Cl^, yep, and Cd^/Cl^ with changes in the vortex lattice. 
However, there is a very significant variation in SCm/^CL (fig- 12). Two different 
types of vortex patterns were employed to produce these variations. The first type used 
uniform values of Nc at each row of horseshoe vortices on the wing-body and on the 

tail. These Nc values were used in combination with three values of Ng. The results 
with uniform distribution of Nc reveal a large variation of BCm/^CL with increasing 
Nc- These results can be shown, by cross -plotting, to be similar to those in figure 7 
because increasing Ng for a given value of Nc has little effect on 9Cm/9CL but 
increasing Nc caused noticeable changes between asymptotic levels of aCm/^CL for 

all values of Ng considered, especially ^ the smaller values of Nc. The second type 

of vortex pattern used uniform values of Nc on the outboard wing panel and outboard 
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tail panel and then used an increased density of elemental panels on the inboard portion 
of the planform. The increased density is illustrated in the input data for sample case 2 
The purpose of these additional inboard elemental panels was to make their chords more 
uniform. This type of vortex pattern virtually eliminated the variation of 

with N(.. These computed results agree with unpublished experimental data for this 
configuration to within O.Olx/cref indicate that good results can be obtained for 
complex planforms with large changes in chord by arranging the pattern of elemental 

panels so that the largest panel chords are no more than two to three times the smallest 
panel chords. 


SAMPLE CASES 


Sample cases have been prepared to illustrate most of the program options avail- 
able. Sketches of the sample cases along with corresponding input data and output data 
listings are provided in appendix C. The sample cases are as follows: 


Sample case Configuration 
1 70 


Description 

Fixed sweep wii^ with dihedral and twist and 
camber 


Page 

46 


2 

3 

4 

5 

6 
7 


13 

113 

no 

15 

215 

315 


Wing-body-tail combination with variable 48 

Wing-body-tail combination with variable 48 

and tail incidence of - 10° 

Wing-body-tail combination with variable sweep 48 

of wing outer panel 

Cropped double-delta wing with variable Ng and 50 
twist and camber to illustrate drag polar option 

Cropped double-delta wing to illustrate C^p 50 

computation 

Cropped double -delta wing to illustrate Cl^ 50 

and Cniq computation 


CONCLUDING REMARKS 


A FORTRAN computer program for estimating the aerodynamic characteristics of 
lifting surfaces in subsonic compressible flow has been described along with the input and 
output variables. Also, a detailed description of the program organization and programed 
equations has been given. The program has been used to compute the aerodynamic 
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characteristics for several configurations that were selected to show the range of plan- 
forms to which the program may be applied. In addition, results from parametric studies 
of the effects of vortex-lattice arrangement on some of the computed aerodynamic char- 
acteristics are presented. From these results, the following recommendations are pro- 
vided as guidance in determining the number of spanwise rows of horseshoe vortices and 
the number of horseshoe vortices chordwise in each row to use to represent a simple 
wing planform or to represent a more complex planform such as a wing-body-tail 
combination: 

1. For simple planforms, (a) use at least 20 spanwise rows and four horseshoe 

vortices chordwise for good values of Cj_,^ ycp> ^D,i/^L^> 

(b) use a vortex-lattice arrangement which gives similar answers for and Cq ^ 

inasmuch as a desirable vortex-lattice arrangement for good values of Cp ji, Cx, and 
Cs is difficult to determine because it is very dependent on the planform. 

2. For a rolling planform, use a large number of spanwise rows and at least two 
horseshoe vortices chordwise. 

3. For a pitching planform, use eight to 10 spanwise rows and eight or more horse- 
shoe vortices chordwise. 

4. For wing -body-tail combinations, use at least 10 to 15 spanwise rows and vary 
the number of horseshoe vortices chordwise so that the local panel chords differ by no 
more than a factor of 2 to 3 from the smallest to the largest. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., October 28, 1970. 
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INPUT DATA 
GROUP ONE 

The input data required for the reference planform is described in the order that it 
is called for by the computer program. All coordinates and slopes should be given for 
the left half of the wing planform. The axis system used is given in figure 1. The y = 0 
intercept coincides with the root chord and is positive pointing along the right wing. 
Although the x = 0 Intercept usually coincides with the intersection of the leading edge 
at the root chord, it may lie anywhere along the root chord; X is positive pointing into 
the wind. All the cards use a format of (8F10.6) for group one data. 

Data for the first card are to be supplied in the following order: 

PLAN Number of planforms for the configuration; use 1 or 2 

TOTAL Number of sets of group two data specified for the configuration 

CREF Reference chord of the configuration 

This chord is used only to nondimensionalize the pitching- 
moment terms and must be greater than zero. 

SREF Reference area of the configuration 

This area is used only to nondimensionalize the computed output 
data such as lift and pitching moment and must be greater than 
zero. 

The data required to define each planform are then provided by a set of cards. The 
initial card in this set is composed of the following data: 

AAN (IT) Number of line segments used to define left half of a wing 

planform (does not include plane of symmetry) 

A maximum of 24 line segments may be used. 

XS (IT) X location of the pivot; use 0 on a fixed wing 

The axis system used is given in figure 1. 

YS (IT) y location of the pivot; use 0 on a fixed wing 

RTCDHT (IT) Vertical distance of particular planform being read in with 

respect to the wing root chord height; use 0 for a wing 
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The rest of this set of data requires one card for each line segment used to define 
the basic planform (variable AAN (IT)). All data described below are required on all 
except the last card of this set; the last card uses only the first two variables in the fol- 
lowing list: 


XREG (I, IT) 


YREG (I, IT) 
DIH (I, IT) 


X location of ith breakpoint 

The first breakpoint is located at the intersection of the 
left wing leading edge with the root chord. They are num- 
bered in increasing order for each intersection of lines in 
a counterclockwise direction. 

y location of ith breakpoint 

Dihedral angle (degrees) in Y-Z plane of line from break- 
point i to i + 1; positive upward 

Along a streamwise line, the dihedral angle is not defined; 
use 0 for these lines. 


AMCD The move code 

This number indicates whether the line segment i is on 
the movable panel of a variable-sweep wing. Use 1 for a 
line which is fixed or 2 for a Une which is movable. 


GROUP TWO 

Three sections of data may be used for group two data. The first section must 
always be included; it is a single card which describes the details of the particular con- 
figuration for which the loading is desired. This card requires a format of (8F5.1, 
F10.4, F5.1, F10.4). The second section is required when the number of horseshoe vor- 
tices used in each chordwise row is not the same; it consists of two or more cards. The 
third section is used when the wing has a twist and/or camber distribution and may con- 
sist of up to 15 cards, depending on the number of horseshoe vortices. The cards in the 
second and third sections use a format of (8F10.4). 

Section one data are to be supplied in the following order: 

CONFIG An arbitrary configuration number which may include up to four 

digits 

sew The number of chordwise horseshoe vortices to be used to 

represent the wing;^ maximum value of 20 may be used 
If set to 0, then a table of the number of chordwise horseshoe 
vortices from tip to root must be provided as TBLSCW (I). 
This sew = 0 option can be used only on wings without 
dihedral and for coplanar wing-tail combinations. 
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VIC 


MACH 


CEDES 


PTEST 


QTEST 


TWIST (1) 


SA(1) 


The nominal number of spanwise rows at which chordwise 
horseshoe vortices will be located 

The variable VIC must not cause more than 50 spanwise rows 
to be used by the program to describe the wing. In addition, 
the product of SSW and SCW cannot exceed 120. If SCW is 0, 
then the sum of the values in TBLSCW (I) cannot exceed 120. 
The use of the variable VIC is discussed in detail in Part I, 
Section 3 of the Program Description. 

Mach number 

Use a value other than 0 only if the Prandtl-Glauert compres- 
sibility correction factor ^ = ^1 - Moo^ is to be applied. It 
should be less than the critical Mach number. 

Desired lift coefficient 

The number specified here is used to obtain the span load 
distribution at a particular lift coefficient. If this answer is 
not required, use 1 for this quantity. If a drag polar for Cj^ 
values from -0.1 to 1 is desired, use 11 for this quantity. 

C;„ indicator 

t-p 

If the damping-in-roll parameter is desired, use 1 for this 
quantity. Except for the incremental pressure coefficients 
and C^p, all other aerodynamic data will be omitted. Use 0 
if C^p is not desired. 

CLq and C^^ indicator 

If these stability derivatives are desired, use 1 for this quan- 
tity. Except for ACp, Cj^, and Cmq, all other aerody- 
namic data will be omitted. It should be noted that both 
PTEST and QTEST cannot be set equal to 1 for a particular 
configuration. Use 0 if Cl^ and Cmq are not desired. 

Twist code for first planform 

If this planform has no twist and/or camber, use a value of 0. 
When this planform has twist and/or camber, use a value of 1 
for this code and provide data for section three. 

Variable-sweep angle for the first planform 

Specify leading-edge sweep angle (degrees) for the first mov- 
able line adjacent to the fixed portion of the planform. For a 
fixed planform, this quantity may be omitted. 
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TWIST (2) Twist code for the second planform 

SA (2) Variable -sweep angle for the second planform 

Section two data are required if SCW is 0. Data for the first variable go on the 
first card and data for the second variable go on the second and following cards. The 
data to be supplied are 

STA Total number of spanwise rows of horseshoe vortices per 

semispan 

This variable sets the number of values of TBLSCW (I) to be 
read in. 

TBLSCW (I) Number of horseshoe vortices in each row starting at the row 

near the tip of the first planform and proceeding to the row 
near the root 

If a second planform has been specified, the table of chord- 
wise rows concludes with number of horseshoe vortices in 
each row of the second planform. For an example, see 
sample case 2. 

Section three data are described as follows: If the configuration has no twist 
and/or camber, the local angles of attack are not specified since the program will set 
them equal to 0. If the configuration consists of two planforms, local angles of attack 
may be specified for both or only one of the two planforms. The twist code describes 
the input to the computer. 

ALP (NV) Local angles of attack in radians 

These are the values at the control point for each horseshoe 
vortex on the wing when the root-chord angle of attack is 0°. 
These data will usually require several cards. For the first 
value on the first card, use the local angle of attack for the 
horseshoe vortex nearest the first planform leading edge at 
the tip; for the second value, use the angle of attack for the 
horseshoe vortex immediately behind in a chordwise direc- 
tion. Continue with the rest of the chordwise row of horse- 
shoe vortices at the tip; then continue inboard at the next 
chordwise row in the same manner to the root until local 
angles of attack for all the control points have been specified. 
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OUTPUT DATA 

The printed results of this computer program appear in two sections: geometry 
data and aerodynamic data. 


GEOMETRY DATA 


The geometry data are described in the order that they are found on the printout. 
The first group of data describes the basic planform; It states the numbers of lines used 
to describe the planform, root chord height, and pivot position and then lists the break- 
points, sweep and dihedral angles, and move codes. These data are a listing of the input 
data except for the sweep angle which is computed from the input data. 

The second group of data describes the particular plarJorm for which the aerody- 
namic data are being computed. Included are the configuration number, the sweep posi- 
tion, a listing of the breakpoints of the wing planform (x, y, and z), the sweep and 
dihedral angles, and the move codes. These data are listed primarily for variable -sweep 
wings to provide a definition of the planform where the outer panel sweep is different 
from that of the reference planform. 

The third group of data presents a detailed description of the horseshoe vortices 
used to represent the planform. These data are listed in nine columns with each line 
describing one elemental panel of the wing in the same order that the twist and/or cam- 
ber angles of attack are provided. (See ALP (NV) in appendix A.) The following items 
of data are presented for each elemental panel; 


X C/4 
X 3C/4 


Y 

Z 

S 

C/4 SWEEP 
ANGLE 


X location of quarter-chord at horseshoe vortex 
midspan 

x location of three-quarter-chord at horseshoe vortex 
midspan 

This is the x location of the control point, 
y location of horseshoe vortex midspan 
z location of horseshoe vortex midspan 
Semiwidth of horseshoe vortex 

Sweep angle of quarter-chord 
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DIHEDRAL 

ANGLE 


Dihedral angle of elemental panel 


LOCAL ALPHA 
IN RADIANS 

DELTA CP AT 
DESIRED CL = 

The fourth group 

REF. CHORD 

C AVERAGE 

TRUE AREA 

REF. AREA 
B/2 

REF. AR 
TRUE AR 
MACH NUMBER 


Local angle of attack at control point (X 3C/4,Y,Z) 

ACp for each elemental panel when wing lift is Cl (j 

of data presents the following geometric data: 

Reference chord of wing 

Average chord (true planform area divided by true span) 

True area computed from planform listed in second 
group of geometry data 

Reference area 

True semispan of planform listed in second group of 
geometry data 

Reference aspect ratio computed from reference plan- 
form area and true span 

True aspect ratio computed from true planform area 
and true span 

Mach number 


AERODYNAMIC DATA 


The aerodynamic data are described in the order that they are found on the print- 
out. Note that Cl^ ^L^Twisf Cxao, Cp^i/CL^, and CL^d are based on 

the specified reference dimensions. 

DESIRED CL Desired lift coefficient specified in input data for 

complete configuration 

COMPUTED ALPHA Cl d/^L(j> angle of attack where desired Lift coeffi- 

cient is developed 

CL(WB) That portion of desired lift coefficient developed by 

the planform with the maximum span when two 
planforms are specified 

When one planform is specified, this is the desired 
lift coefficient. 
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CDI AT CL(WB) 


CDI/(CL(WB)**2) 

1/(PI*AR) 

CL ALPHA 
CL(TWIST) 

ALPHA AT CL = 0 

Y CP 

CM/CL 

CMC 


Induced drag coefficient for lift coefficient in previous 
item 

When two planforms are specified, this is the 
induced drag coefficient of only the planform with 
the maximum span. This result is based on the 
far-field solution (see Part HI, Section 3). 

Induced drag parameter computed from the two pre- 
vious items 

Induced drag parameter for an elliptic load distribu- 
tion based on reference aspect ratio 

^ift-curve slope per radian 

'^ift-curve slope per degree 

Lift coefficient due to twist and/or camber at zero 
angle of attack 

Angle of attack at zero lift in degrees 

Nonzero only when twist and/or camber is specified 

Spanwise distance in fraction of semispan from root 
chord to center of pressure on left wing panel 

Longitudinal stability parameter based on a moment 
center about Y-axis 

Pitching-moment coefficient at Cj^ = 0 


At each chordwise row of horseshoe vortices the following data are presented: 


2y/B Location of midpoint of each chordwise row of horseshoe vortices 

in fraction of semispan locations are listed sequentially from near 
left wing tip toward root 


The next two columns of data describe the additional (or angle of attack) wing loading at a 
lift coefficient of 1 (based on the total lift achieved and the true wing area). 


SL COEF 

CL RATIO 

C RATIO 

LOAD DUE 
TO TWIST 


Span-load coefficient, c^c/ C lc^y 

Ratio of local lift to total lift, c^^Cl 

Ratio of local chord to average chord, c /cav 

Distribution of span-load coefficient due to twist and 
camber at 0° angle of attack 
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ADD. LOAD AT 
CL = 


BASIC LOAD AT 
CL = 0 


SPAN LOAD AT 
DESIRED CL 

SL COEF FROM 
CHORD BD VOR 


Distribution of additional span-load coefficient required 
to produce zero lift when combined with lift due to 
twist and camber 

This distribution is computed at Cj_^ 

Basic span-load-coefficient distribution at zero lift 
coefficient 

These data are the sum of the previous two columns 
of data. 

Distribution of combination of basic span load and 
additional span-load coefficients at desired Cl 

Portion of span-load coefficient due to lift along chord- 
wise bound vortices averaged at horseshoe vortex 
midspan 


In addition, at each chordwise row of horseshoe vortices, the following data are 
presented for induced drag, leading-edge thrust, and suction coefficient characteristics 
computed at an angle of attack of 1 rad from a near -field solution for the additional 
loading (see Part HI, Section 3). 


L. E. SWEEP ANGLE 
CDH C/2B 

CT C/2B 

CS C/2B 

CDH 

CT 

CS 


Leading-edge sweep angle in degrees 

Nondimensional section induced -drag- coefficient 
term 

Nondimensional section leading-edge thrust- 
coefficient term 

Nondimensional section leading-edge suction- 
coefficient term 

Contribution to total drag coefficient from each 
spanwise row of horseshoe vortices, 

Cd,ii(2s cos 4 >) 

Contribution to total leading-edge thrust coeffi- 
cient from each spanwise row of horseshoe 
vortices, ct(2s cos 0)/(q^Sref) 

Contribution to total suction coefficient from each 
spanwise row of horseshoe vortices, 

Cg(2s cos <^) Aq^Sref) 
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Finally, the total coefficient values are listed. 

CDn/CL**2 Induced-drag parameter computed from near- 

field solution 

Leading-edge thrust coefficient computed at 
1 rad angle of attack 

Qg Leading-edge suction coefficient computed at 

1 rad angle of attack 

THIS CASE IS FINISHED End of output for a particular configuration 


For the case where PTEST is 1, all the foregoing aerodynamic output data are 
omitted and only CLP is printed. 

For the case where QTEST is 1, all the foregoing aerodynamic output data are 
omitted and only CMQ and CLQ are printed. 
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SAMPLE CASES 


Input data, sketches, and output data for the sample cases described on page 34 are 
presented in the following order: 


Sample case 

Configuration 

Item 

Page 

1 

70 

Input data 

46 

1 

70 

Sketch 

47 

2,3,4 

13,113,110 

Input data 

48 

2,3,4 

13,113,110 

Sketch 

49 

5,6,7 

15,215,315 

Input data 

50 

5,6,7 

15,215,315 

Sketch 

51 

1 

70 

Output data 

52 

2 

13 

Output data 

59 

3 

113 

Output data 

67 

4 

110 

Output data 

74 

5 

15 

Output data 

80 

6 

215 

Output data 

86 

7 

315 

Output data 

89 


These sample cases reflect the fact that the central processing time for a case is 
generally proportional to the square of the number of horseshoe vortices used to repre- 
sent the left half of a planform. Some typical times for the sample cases with a Control 
Data 6600 computer system are as follows: 


Sample case 

Number of horseshoe vortices 

Time, sec 

1 

100 

62.6 

2 

89 

28.7 

3 

89 

28.7 

4 

52 

7.4 

5 

61 

12.1 

6 

57 

9.0 

7 

96 

34.8 
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STATIC LONGITUDINAL AERODYNAMIC COEFFICIENTS ARE COMPUTED 
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THIS CASE IS FINISHED 


CURVE 5 IS SWEPT 24.73400 DEGREES QN FLANFORM 

CURVE 1 IS SWEPT 0-00000 DEGREES ON PLANFCRM 
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TABtE OF HORSESHOE VORTICES IN EACH CHORDHISE ROW (FROR TIP TO ROOT BEGINNING WITH FIRST PLANFORMI 
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STATIC LONGITUDINAL AEKOOYNAMIC COeFFICIEMS ARE COMPUTED 
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FORTRAN PROGRAM USTING 

This program was written in FORTRAN IV language, version 2.3, for the Control 
Data series 6000 computer systems with the SCOPE 3.0 operating system and library 
tape. Minor modifications may be required prior to use with other computers. The 
program requires 65,000g words of storage on the Control Data 6600 computer system 
and consists of the main program, three overlays, and four subroutines. Each program 
or subroutine is identified in columns 73 to 76 by a 4-character identification. In addi- 
tion, each of these parts is sequenced with a 4-digit number in columns 77 to 80. The 
following table is an index to the program listing: 


Program or subroutine 

Identification 

Page 

WINGTL 

MAIN 

95 

INFSUB 

INFS 

96 

GEOMTRY 

GEOM 

97 

MATXSOL 

MATX 

106 

AERODYN 

AERO 

108 

CDICLS 

CDIC 

116 

MATINV 

MINV 

118 

FTLUP 

TLUP 

120 
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□VERLAy IMINGTL.O.OI 




MAIN 

10 

PROGRAM UINGTU INP UT .OUTPU T , T AI»E5= MPUf , T APE6=0UTPUT ) 


MAIN 

20 

COMMON/ ALL/ BOT r M. SET A> P TE ST >QTES T .TBLSCMIS0)taU20l, 

PN1120! , 

MAIN 

30 

1 PVl 1201 ,AL?l 12 01 ,S 11201 ,?SU 120 J, Prill 1201 

,2Hl 50! 

MAIN 

40 

COMMON/TOTHRE/ CIR ( 120,2 ItSECTRST 1501 




MAIN 

50 

COMMON/ ONETriRE/ TWIST 12 1 .CREF.SREF, CAVE, CLDESfSTRUE 

,AR 

rARTRUE, 

MAIN 

60 

1 RTCOHTI2J, CONFIG, NSSHSVI2I ,MSV 1 2 1 , KBOT , PLAN, 

[PLAN, 

MACH 

MAIN 

70 

2 ,SSUMAI50I 




MAIN 

80 

COMHON/MAINONE/ICOOEOF.TOT AL,AANI2i ,XSI2I ,YS(2I 

,KFCTS (2 1 

MAIN 

90 

1 ,XREGI25,2I . VREGI25,2I ,AREGI25. 21 ,01 HI 2 5, 

21, 

m:D(25, 21 

MAIN 

100 

2 ,XX I25,21,VV (25, 21, AS 1 2 5 , 21 , T TWDI 25 

f2» 

,MMCOI25,2) 

MAIN 

110 

3 ,ANI2I,ZZ 125,21 




MAIN 

120 

7 FORMATllril//lOX,I6,*riORS£SHOE VORTICES LAIOOUT, 

THI S 

IS MORE THAN 

MAIN 

130 

ITHE 120 MAXIMUM. THIS CONFIGURATION IS ABORTED. 

♦! 



MAIN 

140 

8 FORMAT IlHl// lOX, 16 ♦ ROWS OF HORSESHOE VORTICES 

LAIDOUT. THIS 

IMAIN 

150 

IS MORE THAN THE 50 MAXIMUM. THIS CONF 15URATI0N 

IS 

ABORTED.* 1 

MAIN 

160 

9 FORMAT IlHl // lOX, *PLANFORM* 16 ♦ HAS* 16 




MAIN 

170 

1 ♦ BREAKPOINTS. THE MAXIMUM OIMENSIONEO IS 25. 

THE 

CONFIGURATION 

IMAIN 

ISO 

2S ABORTED. ♦! 




MAIN 

19 0 





MAIN 

200 

VORTEX LATTICE AERODYNAMIC COMPJTATION 




MAIN 

210 

NASA-LRC PROGRAM NO. A279^ 




MAIN 

220 





MAIN 

230 





MAIN 

240 





MAIN 

250 

ICODEOF=TOTAL=0 




MAIN 

260 

WINGTL^6LMINGTL 




MAIN 

270 

RECALL=6HRECALL 




MAIN 

280 

1 CALL OVERLAYIWINGTL, 1,0, RECALL! 




MAIN 

290 

IFIICODEOF.GT.O) GO TO 99 




MAIN 

300 

1FIM.GT.120! GO TO 2 




MAIN 

310 

NS^ = NSSWSVm NSSKSVI2I 




MAIN 

32 0 

IF ( NSM*GT.50 ! GD TO 4 




MAIN 

330 

ITSV = 0 




MAIN 

340 

DO 10 IT=1,IPLAN 




MAIN 

350 

IF i AN! ITI.LE.25. 1 GD TO ID 




MAIN 

360 

WRITE (6,91 IT,ANIITi 




MAIN 

37 0 

ITSV = 1 




MAIN 

38 0 

10 CONTINUE 




MAIN 

39 0 

IF (ITSV. GT. 01 GO TO 5 




MAIN 

400 

GO TO 3 




MAIN 

410 

4 WRITE 16,81 NSW 




MAIN 

42 0 

GO TO 5 




MAIN 

430 

2 WRITEI6,7T M 




MAIN 

440 

GO TO 5 




MAIN 

450 

3 CALL OVERLAY! WINGTL ,2 ,0, RECALL 1 




MAIN 

46 0 

CALL OVERLAYIWINGTL, 3,0, RECALL! 




MAIN 

470 

5 TOTAL=TOTAL-l. 




MAIN 

48 0 

GO TO 1 




MAIN 

490 

99 STOP 




MAIN 

500 

END 




MAIN 

510 
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SUBROLTINE INFSUB ( 00 T * Ft I ♦ F V 1 1 FW II 

CGMMCN/ INS UB2 3/PS IIfAPHlIt>XX,YYYiZ22fSNNfTQLFNC 

FC =CCS(PSin 

FS =SIN(PSm 

FT =FS/FC 

FPC=CCSIAPHII I 

FPS = SIN(APhI I ) 

FPT=FPS/FPC 
FI =>XX<-SNN*FT*FPC 
F2 =YYY'fSNN*FPC 
F3 =ZZZ^SNN*FPS 
F^ =XXX-SNN*FT*FPC 
F5 *VYY-SNN*FPC 
F6 =ZZZ-SNN*FPS 

FFA= ( X XX* *2 ♦( VTV *F PS )♦*2+F pc* *2*( ( YYY*FT) **24 (ZZZ/FC )**2-2.* 

1XXX*YYY*FT)-2.*ZZZ*FPC*I VYV*FPS+XXX*FT*FPSI ) 
FFB=(F1*FHF2*F2^F3*F3)**.5 
FFC=(FA*F4*F5*F5*F6*F6)*«.5 
FFD=F5*F5+F6*F6 
FFE=F2*F2^F3*F3 

FFF=( F1*FPC*FT+F2*FPC + F3*FPS ) /FF0 - ( F4*F PC*F T4 F5*F PC4-F6*FPS ) / 

IFFC 


TKG TCL6KAXCE SET AT THIS PCINT IN THE PROGRAH MAY NEEC TO BE 
CHANGED FOR CCHPUTERS OTHER THAN THE COC 6000 SERIES 


IF(ABS(FFA).LT.(BCT*L5.E-5)**2) GO TO 262 
FUUNE=(ZZZ*FPC-YYY*FPS)*FFF/FFA 
FVONE=IXXX*FPS-ZZZ*FT*FPC)*FFF/FFA 
Fk*CNE= I YYY*F T-XXX )*FFF/FFA*FPC 
GO TO 265 

262 FUGNE=FVnNE=FWCNE=0. 

265 IF (ABS(FFD) .I.T.TOLBNC) GC TO 263 

FYTWO= F6* (U-F4/FfC)/f FC 
FWTW0=-F5* (1.-F4/FFC ) /FFC 
GO TO 266 

263 FVTW0=FkTWC=0. 

266 IF(ABS(FFE).LT.T01RNC) GC TO 264 

FVTHRE=-F3*( l.-Fl/FFB)/FFE 
FWTHRE=F2* ( 1 .-F 1/ F F0 I /FF E 

GC TC 267 

264 FVTHRE^FWTHRE-C, 

267 FUI=FUCNE 
FVI=FVCNE4FVTWC+F VTHRE 
FWI = FWCNe + FV^ThO+FWTHRE 
RETURN 

END 


INFS 10 
INFS 20 
INFS 30 
INFS 40 
INFS 50 
INFS 60 
INFS 70 
INFS 80 
INFS 90 
INFS 100 
INFS 110 
INFS 120 
INFS 130 
INFS 140 
INFS 150 
INFS 160 
INFS 170 
INFS 180 
INFS 190 
INFS 200 
INFS 210 
INFS 220 
INFS 230 
INFS 240 
INFS 250 
INFS 260 
INFS 270 
INFS 280 
INFS 290 
INFS 300 
INFS 310 
INFS 320 
INFS 330 
INFS 340 
INFS 350 
INFS 360 
INFS 370 
INFS 380 
INFS 390 
INFS 400 
INFS 410 
INFS 420 
INFS 430 
INFS 440 
INFS 450 
INFS 460 
INFS 470 
INFS 480 
INFS 490 
INFS 500 
INFS 510 
INFS 520 
INFS 530 
INFS 540 
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Overlay jwiNGTLt I »o) 

'coB»ON/H«!NoiE/ICOOEOF.tOT«I..A«N(Jl.XSI2II.YSi;2I.KFCISI|l 

3 ,AN(2)fZZ (25*2) 

1 F0RMAT*nHl// 63X,*6E0METRY DATA* ) 

2 FORMAT </// 45X .AlO, *«P^ERENCE P^ANFORM HAS 

I 12X -ROOT CHORD HEIGHT »* Pl|.5 » AX 

7 PIVOT POSITION* AX *X(S) ** F12.5*5X *Y(S» »* 

3 *0REAK POINTS FOR THE REFERENCE PLANFORH 

3 FORMAT I8F10.A) 

5 Format *aHi*// atx « *configuRation no.* »f8.o / > 

6 FoRMAT(2F12.5,2E12.5,F12.5) 

7 FORMAT! 


GEOH 

10 

GEOH 

20 

GEOH 

30 

GEOH 

40 

GEOH 

50 

GEOH 

60 

GEOH 

70 

GEOH 

80 

GEOH 

90 

GEOH 

100 

GEOH 

110 

GEOH 

12Q 

GEOH 

X3Q 

GEOH 

140 

GEOH 

150 

GEOH 

160 

GEOH 

170 


13 * CURVES* // 

*VARIABLE SWEEPGEOM 
F12.5 //A6X» 

* / ) 


180 


6E0H IRO 
GEOV 200 
GEOH 210 
GEOM 220 
GEOH 230 
GEOH 2A0 

, „™.. . OFSESHOE FOBj;iCES^O^ LEFT HALF OF TOE aO|OF |E. 

1In6/36X,IA»10H CH0RDWISE,21X,IA,9H SPANWISE/^ *SmEEP* 7X *DIHEDRAGE0H 270 

8 format !22X *P0INT* 6X *X* UX *Y* UX *Z lOX SHEtK 
IL* AX *MOVE* / 68X *ANGLE* 8X *AnGLE 6X CODE / 

9 format (20X, IS, 3F12.5.2FU.5, 16) , DEGREES Oh PLANFORGEOH 300 

10 EORMAT ( / AOX, ‘CURVE* 13 * IS SWEPT* 1-12.3 uco c 

n*F 0 RMAT(lHl///AlX •rct^vARIABlE^SWEEp'^OL^ =*'gEOH 330 

,E^FoR-AT iS^r^toE^F^E .A, OEO- 3A. 

13 FORMAT !8F5.1,F10.A,F5.1,F10.A) 

1? FSR5«'a5i'///”^*E«"0»'-’'’"““*" PROCESS PTEST F5.1 

15 ;5EF":or:.J"E. nr.AiSfE. ,x -code. / 

!S fr.:^x}xXERx!5A;T'HrS2fAK:fe;S”[Sc5;irs'pSrw5ll A^.^..S.3X.E0HFOEC. ^ 

IAS BEEN adjusted qE yORTICES IN EACH CHORDWISE ROW* ) GEOH A50 

IS Fsr^xllx.^Jl;E"oF^ SI 5 V ? -"'iis: :ss 

EA'ioii^rxxr3s;.ii';isisiHiE >si;iis p 

1OnFIGURATION*//S0X*PLANFoRM total SPANWlbE 

25 Format (szx, ia , lox , i3 » iix , ia ) 


GEOW 3A0 
GEOM 350 
GEOM 360 
GEOM 370 
GEOM 380 
GEOH 390 
GEOH 400 
)GEOH 410 
GEOH 420 


PART ONE - geometry COMPUTATION 

SECTION ONE - INPUT OF REFERENCE WING POSITION 


RTCOHT ( 1 ) sRTCOHT ( 2 ) =0 . 
YTOL = UE-10 


GEOH 500 
GEOH 510 
GEOH 520 
GEOH 530 
GEOH 540 
GEOH 550 
GEOH 560 
GEOH 570 
GEOH 580 
GEOH 590 
GEOH 600 
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c 

c 

c 


c 

c 


AZV = l.E+13 

PIT « 1.5707963 

rad « 57.29578 

IF ( TOTAL . GT.O. ) GC TO 80 


SET PLAN EQUAL TO 1. 
VARIABLE SWEEP WING 
SET PLAN ECUAL TO 2. 


FOR A WING ALCNE COMPUTAICN - EVEN FOR A 
FOR A WING - TAIL C0M8INATICN 


SET TOTAL EQUAL TC THE MjPeER CF SETS 
OF GROUP TWO DATA PRCVICEC 


READ <5, 3) PLAN, TOTAL, CREF.SREF 
IF «ECF,5| 1CC6,AC 
40 IPLAN =PLAN 


^QUAL TO TEE KA)(IHLP NUPBER OF CURVES PEQUIREC TO 
define the PLANFORM PERIFETEP CF THE (IT) PLANFCRP. 

SET RTCDHTIIT) EQUAL TO THE ROCT CHORD HEIGHT CF TEE L IFTiKr 

iitkwhose peripeter points are being rL; in' wl^H 

RESPECT TO THE WING ROOT CHCRO HEIGHT 

WRITE (6,1) 

DC 58 IT = 1, IPLAN 

READ (5,3) AAN( IT J , )<S { I T ) , VS < I T ) , RTCCFT C IT J 
N = AANUT) 

N1 = N ♦ 1 

MAK = 0 

IF (IPLAN. EQ.l) PRTCCN = lOH 

.and. IT.EQ.l ) PRTCON = lOH 

JD » 1 ^ «AND. IT.EO.2 ) PRTCON = lOh 

WRITE ( 6,21 PRTCON, N,RTCCHT(IT),XS(IT),YS(m 
WRITE(6,17) 

DC 59 1=1, N1 

Ki:m’= ciH.i.in, APCD 

IF { I .EO. 1 ) GO TO 59 

IF ( WAK.NE.O .or. HCC ( I-l , IT) .NE.2 } GO TO 49 

MAK = I—l 

GC TO 51 

50 YREGII.IT) = VREG(I-1,IT) 

AREG( I-1,IT) = AZY 

ASWP < 90. 

51 J =1-1 


FIRST 

SECOND 


WRITE PLANFCRP PERIRETER PCINTS AND ANGLES 

59 CONTINUE 

KFCTSIIT) = HAK 

58 cSIunIe’^"’ M.XREG«N1.IT),YREG(N1,IT) 


part 1 - SECTION 2 


GEOM 610 
GEOM 620 
GEOM 630 
GEOM 640 
GEOM 650 
GEOM 660 
GEOM 670 
GEOM 680 
GEOM 690 
GEOM 700 
GEOM 710 
GEOM 720 
GEOM 730 
GEOM 740 
GEOM 750 
GEOM 760 
GEOM 770 
GEOM 780 
GEOM 790 
GEOM 800 
GEOM 810 
GEOM 820 
GEOM 830 
GEOM 840 
GEOM 850 
GEOM 860 
GEOM 870 
GEOM 880 
GEOM 890 
GEOM 900 
GEOM 910 
GEOM 920 
GEOM 930 
GEOM 940 
GEOM 950 
GEOM 960 
GEOM 970 
GEOM 980 
GEOM 990 
GEOMIOOO 
GEOMlOlO 
GEOM1020 
GE0M1030 
))GE0M1040 
GE0M1050 
GE0M1060 
GEQM1070 
GE0M1080 
GE0M1090 
GEOMllOO 
GEOMlllO 
GE0MU20 
GE0M113G 
GEOMl 140 
GE0M1150 
GEOM1160 
GE0M1170 
GE0M1180 
GE0M1190 
GEQM1200 
GE0M1210 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


ftE^C GPCUP 2 C^TA ANC CCFPUTE DESIRED WING PCSITICN 


GEOM1220 
GE0M1230 
GE0M1240 

sew MUST NOT BE SET EQUAL TC ZERC CR CNE WHEN THE WING HAS DIHECRALGEOM1250 

GEOMl 260 

SET SA(l)fSA(2i EQUAL TC THE SWEEP ANGLEfIN DEGREES# FCR THE F IRSTGE0M1270 
CUkVEIS) THAT CAN CHANGE ShEEP FOR EACH PLANFCRM GEOM1280 

GEOMl 290 

IF A PARTICULAR VALUE OF CL IS DESIRED AT WHICH TFE LOADINGS ARE GEOM1300 
TO BE COMPUTED, SET CLOES EQUAL TO THIS VALUE GEOM1310 

SET CEDES EQUAL TC 11. FCR A DRAG POLAR AT CL VALUES CF-.l TO 1.0GE0M1320 

GEOM 1330 

IF PTEST IS SET EQUAL TC CNE TFE PROGRAM WILL COMPUTE CLP GE0M13A0 

IF OIEST IS SET EQUAL TC CNE TFE PROGRAM WILL COMPUTE CMQ AND CLQGE0M1350 

DC NOT SET BOTH PTEST ANC CTEST TC ONE FOR A SINGLE CONFIGURATION GE0M1360 

GEOMl 37 0 

CR TWIST(2) EQUAL TO C* FOR A FLAT PLANFCRM AND TC 


SET TWIST(l) 

FOR A PLANFCRM THAT HAS TWIST AND/CR CAMBER 


GEQM1390 
GEGM1400 

80 READ ( 5tl3 ICCNF IG #SCW # V I C# MACH#CLDES# PTEST,QTEST # TWIST ( I ) # SA( I ) 

1ISTI2J #SAC 2) 

WRITE(6#5I CCNFIG 
IF (ECF#5) 10C6.82 

82 IF ( PTEST. NE.O. .AND. CTEST. NE.O. ) GO TO 1CC8 
IF (SCW.EC.O.) GC TC 76 

DC 74 1=1 #50 
74 TBLSewm = sew 
GC TO 78 

76 READ (5#3) STA 

NSTA = STA , 

READ (5# 3 1 (TBLSCWn ) #TBLSCW( I ♦ 1 ) # TBL SCW ( I +2 ) » T®*- ^CW( I *3 ) 

1 


IT 


#TBLSCW( 1 + 4 ) #TBLSCW( 1 + 5) f TBLSCWI 1 + 6 I #TBLSCW( 1 + 7) , 
I = 1#NSTA,8) 

1 flPLAN 


AAM IT) 

N + 1 

XREG( I»IT) 
YREGI I, IT ) 
AREGI 1# IT ) 
ATAN( A(l) ) 


RSAPin = PIT 


GO TC 79 


78 DC 100 
N 

Nl 

DC 83 I=1#N 
XREFII) 

YREFII) 

A (I ) 

RSAR(I) 

IF (A( D.EC.AZY) 

83 CCNTINLE 

XREF(M) = XREG(N1,IT) 

YREF(Nl) = YREG(M,IT) 

IF ( KFCTS(IT) .GT. 0 ) 

K =1 

SA( IT) = RSARC 1) ♦ RAO 
GC TO 77 

79 K = KFCTSI IT) 

77 WRITE (6,10) K,SA(IT)#IT 
SB = SA(IT)/RAO 

IF ( A8S( SB - RSARIK) ).GT. (.1/RAC) ) GO TC 111 
REFERENCE PLANFCRM COORDINATES ARE STORED UNCHANGED FDR WINGS 
WITHOUT CHANGE IN SWEEP 
DC 113 I=1#N 
X( I ) = XREF( I) 

Y( I )=YREF( I) 

IF (RSAR(I) .EC. PIT I GO TC 114 

A(n=TAN(RSAR(I)) 

GC TC 113 


GE0M1420 
GE0M1430 
GE0M1440 
GEOM1450 
GEOM1460 
GEOM1470 
GE0M1480 
GE0M1490 
GEOM1500 
GE0M1510 
GEOMl 520 
GE0M1530 
GEOM1540 
GE0M1550 
GE0M1560 
GE0M1570 
GE0M1580 
GE0M1590 
GE0M1600 
GEOM1610 
GE0M1620 
GEGM1630 
GE0M1640 
GE0M1650 
GE0M1660 
GE0M1670 
GE0M1680 
GE0M1690 
GE0M1700 
GE0M1710 
GEOM1720 
GEOM1730 
GEOM1740 
GE0M1750 
GE0M1760 
GE0M1770 
GE0M1780 
GE0M1790 
GE0M1800 
GE0M1810 
GE0M1820 
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114 A(IJ=:^ZY 
113 SARC I )=RSAft( I ) 

X(N1 | = XREF (M ) 

Y(M I=YREF(N1 ) 

GC TO 103 

CHANGES IN HING SWEEP ARE PACE HERE 

111 IF (MCC(K, IT).NE.2) GC TC 1007 

KA=K-1 

QC 81 1=1, KA 
XI I ) = XREFC I) 

Yl I »=YREF ( I ) 

81 SARI I l=RSAR(I ) 

C DETERMINE LEADING EDGE INTERSECTION BETWEEN FIXED AND VARIABLE 

C SWEEP WING SECTIONS 

SAR(K)=SB 
AIK) * TANISB) 

SAI=SB-PSAR(K ) 

XIK + I ) = XS+(XKEF(K*1 )-XS)*CCS(SAI )+(YPEF(K + l )-YS)^ SINI SAI ) 

Y(K<^I) = YS+ (YREFIK41 )- VS ) *CC S I SA I) - ( XREF ( K ♦! )-XSJ*SIM SAI } 

IF I ABS (SB - SARIK-1) ) .LT. I.l/RAO) > GC TC 84 

Y(K)=X(K<'1)-X(K-1)-AIK)*V(K + 1)+A(K-1)*YCK-1) 

Y(K)=Y(K| /( AIK-D-AIKH 

X(K)= AIK) ♦X(K-l)-ACK-l)4X(K + l)4A(K-n*A(K)^(Y{K+n-Y(K-l) ) 
X(K)=X(K)/(A(K)-A(K-U) 

GC TC 65 

C ELIMINATE EXTRANECLS 8RE/KPC1NTS 

84 X(K) = XREF(K-1 ) 

Y(K) = YREF(K-1 ) 

SAk(K) = SARIK-1) 

85 K=K+1 

C SWEEP THE BREAKPOINTS CN TEE VARIABLE SWEEP PANEL 

C (IT ALSC KEEPS SWEEP ANGLES IN FIRST OR FOURTH CUADRANTS) 

86 K=K4l 

SAK(K-1)=SAI4P5AR(K-1) 

99 IF ( SAR(K-l) .LE. PIT ) GO TO 102 

SAR(K-i ) = SAR( K-1)-3.1^I5S27 
GC TO 99 

102 IF ( SAR(K-i) •GE.(-PIT)) GO TC 106 

SARIK-1 )=SAR( K-1 ) 43.14 15S27 
GC TO 102 

106 lf(I SARIK-I) ).LT..C) GC TC 1C8 

IF I SARIK-1) - PIT ) 90,87,87 

108 IF I SAR(K-n 4 PIT ) 89,89,90 

87 AIK-l)=A2Y 
GC TO 91 

89 AIK-1)=-AZY 
GC TO 91 

90 AIK-1)=TANISARIK-1) ) 

91 KK = MCDCK , IT ) 

GC TO 193*92) ,KK 

92 Y(K) = YS4l y REF IK )- YS )*COS (SA I )-( XREF (K )-XS) M SA I) 

X(K)=XS4l XREF IK )- X S)*COS ISA I )4(YREF IK )-YS)^SI MSA I) 

GC TO 86 

C DETERMNE THE TRAILING ECGE INTERSECTION 

C BETWEEN FIXED AND VARIABLE SWEEP WING SECTIONS 

93 IF (ABS (RSAR(K)-SARlK-in .LT. (.1/RAO) ) GC TO 96 

V(K) = XREF (K41 )-X(K-l)-A (K)*VREF(K41 )4A(K-1)*Y(K-1 ) 
Y(K)=Y{K)/(A(K-1)-A(K)) 

X(K) = A(K)*X(K-1 )-A(K-l )*XREF(K41)4A(K-1)*A(K) ♦(YREF(K4 1)-Y(K-1) ) 
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XCKJ=X{K)/<A(K|-ACK-1)) 

GC TO S7 

96 X(K) = XREF(K + n 
Y(K)=YREF<K+1) 

C STORE REFERENCE PLANFCRM CCCRDINATES CN INBOARD FIXED TRAILING 

C EDGE 

OG 9Q I=K.M 

xm=XREF( I) 

Y(n=YREFl I) 

98 SARI I-1)=KSAR(I-1 I 
103 DC 101 I=lfN 

XXI IrIT) = xm 
YY(IflT) = YIII 
MKCDIIfIT)= MCC( I, ITI 
TTWOI I » IT ) = CIH( I f IT ) 

101 AS ( I,IT) = A( I) 

XX(NlflT) - XINI) 

YYlNltfT) = Y(Nl) 

AN(IT) = A AN! IT I 
100 CONTINUE 

LINE UP BREAKPOINTS ARCNG FLANFCRPS 


299 BGTSVI 1) = BCTSV(2)=0. 

WRITE (6,16) 

DO 180 IT=lfIPLAN 
NIT = AM IT) +1 

00 178 ITT=l*IPLAN 

IF ( ITT.EO.IT) GC TC 178 
NITT=AN( ITTI+1 
DC 176 I=lfNlTT 
JPSV=0 

DC 166 JP=ltNIT 

IFIYYUPflT) .EO.YYI If ITT) > GO TC 176 
166 CONTINUE 

DO 17C JP=lfNIT 

IF lYYUPf m.LT.YYI If ITT) ) GC TO 168 
170 CONTINUE 
GC TC 176 
168 JPSV = JP 

IND = NIT -UPSV -n 
DC 172 JP=lfIND 
K2 * MT -JP +2 
K1 = MT -JP 
XX(K2»m = XXlKlflT) 

YYIK2f IT) = YYIKI fIT) 

MKCD(K2fIT)= PPCDIKl.IT) 

AStK2tIT) = AS(KlflT) 

172 TTWDIK2tIT)=TTW0(Kl,IT) 

YY(JPSVflT) = YYdf ITT) 

ASIJPSVfIT) = AS(JPSV-lflT) 

TTWDI JPSV, IT)= TTWCI JPSV-lfIT) 

XX(JPSV,IT) = (YYUPSV.IT) - YY UPSV-1 , IT ) ) 

1 + XXI JPSY-lflT) 

MKCOI JPSVf IT ) = MPCCI JPS\<-1,IT) 

AM IT) = AM IT) ♦ 1. 

nit = MT ♦ 1 
176 CONTINUE 
178 CONTINUE 
C 


♦ AS(JPSV-lflT) 
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SEQUENCE WING COORCINATES FPCM TIP TC ROOT 

Nl = AN(IT)* !• 

OG 203 1=1 , Nl 
203 0< I J = YY( I, IT) 

OC 208 J=1,N1 
HIGH = 1. 

OC 205 1=1 , Nl 

IF {{ C(n-HIGH)*GE.O. ) GO TC 205 

HIGH = on ) 

Ih = I 

205 CONTINUE 

IF U.NE.l) GO TO 2C6 
BCTSVdT) = HIGH 
KFX< IT) = Ih 

206 0 ( IH) = 1. 

SFYUfIT) = hIGH 

IF UF.GT*KFX(IT) ) GO 1C 2C9 
IYL( J, IT) = 1 
IYT( J, IT) = 0 
GO TO 208 
2C9 lYLUflT) = 0 
lYTUtlT) = I 
2C8 CONTINUE 
180 CONTINUE 

SELECT MAXIMUM 8/2 AS TF£ WING SEMISPAN 
KBOT = 1 

IF (BCTSVI l).GE,BOTSV(2>) KBCT = 2 
BCT = eCTSV(KBCT) 

COMPUTE NOMINAL HORSEShCE VERTEX WIDTH ALONG WING SURFACE 


TSPAN = 0 

ISAVE = KFX(KeCT) - 1 
I = KFX(KECT) - 2 


216 

IF lI.EC.O) 

GC 

TO 

217 

217 

IF (TTWCd ,K60T) .EC.TTWDI I SAVE, KBCT) ) 
CIWD = COS( ATANITTWCI ISAVE, KBCT) ) 

GC 

) 

TO 

218 


TLGTh = (YVdSAVE + l.KBCT ) - Vy(I*l,K8CTI 
TSPAN = TSPAN ♦ TLGTH 

) 

/ CTWC 


IF (I.EC.O) 
ISAVE = I 

GC 

TO 

219 


218 I = I -1 

GO TC 216 


219 VI = TSPAN / VIC 
VSTOL = VI / 2 

ELIMINATE PLANFORM BREAKPOINTS WHICH ARE WITHIN (B/2)/2000 UNITS 
LATERALLY 

OC 22C IT » 1 ,IPLAN 
N = AM IT) 

Nl= N ♦ 1 
00 22C J=1,N 

AA = ABS(SPYU.IT) - SPy(j4l,IT) ) 

IF ( AA.EQ.G. .OR, AA .G 1 .AB S (T SP AN/2CC0 . ) ) GC TO 22C 
IF { AA.GT.YTCL) WRITEI6,19) SPYU4l,IT) , SPY(J,IT) 

00 222 1=1 ,N1 

IF ( YY( I , IT ) .NE. SPY! J+1 ,IT ) ) GC TO 222 
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c 

c 

c 


YY( I,m = SPY(J,IT) 

222 CONTINUE 

5FY(J*l,m = SPY<J,IT) 

220 CONTINUE 

COMPUTE I COORDINATES 

00 236 IT^l.IPLAN 

JK * N1 = AN(IT) + 1. 

DC 230 JZ=1,N1 

230 2Z(J2,IT> = RTCOHT(IT) 

JZ = 1 

232 JZ * JZ > 1 

IF UZ.GT.KFXUT)) GC TO 234 

ZZ(JZtlT) = ZZ<JZ-1,IT) 4(VV(JZ,IT) - YY(JZ-1,IT» ) ♦TTWDUZ-l,IT 
GC TO 232 

234 JM » JM-l 

IF ( JM.EU.KFX(IT) » GC TC 236 

ZZ(JHflT) = ZZ(JM4l,m 4<yy(JM,IT)-VY(JM+l.m) ♦TTW0UM,IT) 

GO TO 234 
236 CONTINUE 

WRITE PLANFORM PEfilFETEP PCINTS ACTUALLY USED IN THE CCMPUTATICNS 

WRITE (6»6) 

DC 240 IT =1 , IPLAN 
N = ANIIT) 

N1 » N ♦ 1 

IF (IT.E0.2) WRITE (6,10 
DC 238 KK=1,N 

TOUT = ATAN ( TTWD(KK,IT) )♦ RAD 

ACUT = ATAN( AS(KK, IT) )ARAC 

IF (AS(KK, IT) .EQ.AZY) A0UT=90* 

WRITE (6,9) KK,XX(KK,IT) , YY(KK,IT), 2Z(KK,IT), ACUT, 

1 TOUT tMMCC(KK,IT) 

238 CONTINUE 

WRITE (6,9) N1,XX(M,ITI,YV(N1,IT),ZZ(M,IT) 

240 CONTINUE 

PART ONE - SECTION THREE - LAY CUT YAWED HORSESHOE VORTICES 


STRUE = 0. 

NSSWSV(l) = NSSWSV(2) = HSV(l) = MSV(2) = C 

700 DC 722 IT=1,IFLAN 

N1 = AN(IT) ♦ 1. 

I « 0 

J =1 

YIN = BOTSV(IT) 

ILE = ITE = KFX(IT) 

0 DETERMINE SPANWISE EOROERS CF hCRSESFCE VORTICES 

701 IXL = IXT 0 

I =1*1 

IF(YIN.GE. (SPYC J, IT)+VSTCL) ) GO TO 703 
C BORDER IS WITHIN VORTEX SPACING TOLERANCE (VSTCL) CF BREAKPOINT 

C THEREFORE USE THE NEXT BREAKPOINT INBOARD FOR THE BORDER 

VBOROm * YIN 
GC TO 707 

C USE NOMINAL VERTEX SPACING It DETERMINE THE BORDER 
703 VBORO(I) = SPY(J,IT) 

C COMPUTE SUBSCRIPTS ILE AND ITE TO INDICATE WHICH 

C BREAKPOINTS ARE ACJACENT AND WHETHER THEY ARE CN IKE WING LEADING 
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c 


c 


715 


706 


707 


708 


709 


tOGE CR THE 
IF (J.GE.Nl) 
IF (SPY(J»IT* 


trailing ecge 


.NE. SPY( J^l f IT ) ) 


GO 

GO 


TO 

TC 


706 

706 


IXL 

IXT 

J 

GC TC 

YIN 

IXL 

IXT 

J 

CFHI 

IFHI 

IF ( 

YIN 

IF 

ILE 

ITE 


IXL 
IXT 
J + 


lYLUtlT) 
lYTU, IT) 


715 


= CCS ( 
= ILE - IXL 
J.GE.M ) 

= YIN - VI* 
i.NE.n 

* ILE - 
= ITE 


SPY( J,IT ) 

IXL ♦ lYU Jf IT ) 

IXT ♦ IYT( J, IT) 

J + 1 

ATAN ( TTi^Cl ILE, IT) 


) ) 


IPHI 
CCS ( 


= 1 
/TAN 


TTWOI IFhl.IT) 
GO TO 709 


) ) 


IXL 

IXT 


GG TO 701 

CLMPLTE COCRLINATES 


YC 
Hk 
IKl 

7h{ IMl) 

Phl( IHl) = TTkC( ILE, IT) 

SSWWA(IMl) = ASULEtlT) 
XLE = XXI ILE, IT) ♦ 

XTE == XX(ITE,IT) ♦ 

XLOCAL = ( XLE - XTE 


FCR CHCRDklSE RCk CF HORSESHCE VCRTICES 
VBOROd^l) + VBCRD(I) ) / 2. 

( VBORC(I) - VBCRCII-1))/ 2. 

I - 1 ♦ NSSkSVm 

ZZ(ILE,IT) ♦ ( YC - YYlILEtIT) ) * TTkC<ILE,IT) 


AS( ILE, m ♦ (YQ - 
AS( ITE, IT) * (YQ - 
/ TBLSCwnWl) 


YY( ILE,IT) 
YY(ITEflT) 


COMPOTE WING AREA PROJECTED TO THE X - Y PLANE 

STRUE = STBUE + XLOCAL ♦ TBLSCW(IMl) ♦ (HW ♦ 2.) * 2. 

NSCW = TBLSCk(IMU 

OG 72C JCW=1,NSCW 
AJCW = JCW * 1 

XLEL = XLE - AJCW ♦ XLCCAL 

NTS = JCW ♦ MSV(l) + )^SV(2) 

PMNTS) = XLEL ^ .25 ♦ XLCCAL 

PV(NTS) = XLEL - .75 ♦ XLCCAL 

PSI(NTS) = ((XLE - PN(NTS))*AS( ITEfIT) 

1 IT) ) / (XLE - XTE) ♦ CPHI 

= HW / CPHI 
= YC 


♦ (PN(NTS) 


XTE)*AS( ILE 


720 


S(NTS) 
0(NTS> 
CONTINUE 
MSV( IT) 


= MSV( IT) ♦ NSCW 


TEST TO DETERMINE WHEN WING ROCT 
IF I VBORC( I ) .LT. -C. ) 


(Y = C) 
GO TO 


IS REACHED 
708 


NSSWSV(IT) = I - 1 
722 CONTINUE 

M = MSVd ) ♦ MSV(2 


COMPUTE ASPECT RATIC ANC 


AVERAGE CHORD 


SOT = - BOT 

AR = A. ♦ eOT ♦ ECT / SREF 

ARTRUE = 4. ♦ BCT ♦ BCT / STRUE 

CAVE = STRUE / ( 2. ♦ ECT ) 
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beta = i 1. - MACE* MACE) ♦♦ .5 

NVTWO = 0 

DC 354 IT*1,IPLAN 

iNIVONE • I ♦ (IT-n*MSV(l) 

NVTWC = UMM + MSV(1T> 

IF (ThISKITI .LE. 0. ) GQ TC 350 

READ (5*3) (ALP(NV)fALP(NV+l)»ALP(NV*2) * ALP( NV+3 } » ALP ( NV+4 ) *AL 
I 4^5) *ALP(NV + 6 )*ALP(NV + 7) ,NV=NVONE* NVTWCtB) 

GC TO 354 

350 DC 351 NV = NVCNE , NVThC 

351 ALP(NV) * 0. 

354 CONTINUE 

WRITE (6,24) M 

WRITE (6*25) (IT,MSV(IT)*N5SlrtSV(IT)* IT=1*IPLAN) 

IF C SCW.NE.O* ) WRITE (6*20) SCW 

IF ( SCW. EG. 0. ) WRITE (6*22) ( T BLSCW ( I ) , 1 = 1 * NST A ) 

APPLY PRANOTL-GLAUERT CCRRECTICN 


DC 360 NV = 
PSKNV) 

PN (NV) 

360 PV (NV) 
RETURN 

1006 ICODECF 
WRITE(6*11) 
RETURN 

1007 ICODECF 
WRITE(6,12) 
RETURN 

1008 ICODECF 
WRITE (6*15) 


l,M 

ATANI PSI(NV)/BETA) 
PN(NV) / BETfl 
PV(NV) / BETA 

1 

CONFIG 


2 

K*IT 


3 

PTEST.OTEST 


RETURN 

END 
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OvERLAY(WIN6TL*2»0) 

Program matxsol 

Dimension yy< 2) fPUizj tEVis) trw( 2 » .fvn(i20«i 20) tiPivornao ) * 

1 rNDEX(120«2) 

Common/all/ BOT.M*BETA,PTESTfOTEST.TBLSCH(SO) •OU20> »PNU20> * 
1 PV (120) t ALP ( 120 ) •s(120) «PSI ( 120 ) yPHI < 120 ) .ZH( 50 ) 

COMMON/TOTHRE/ CIR(120.2) ,SECTRST(50) 

C0MM0N/INSUB23/APSI»APHI tXX tYYYtZZ tSNN.TOLC 


Part 2 - compute circulatio)^ terms 


FpI = 12.5663704 


The tolerance set at this point IN THE PROGRA)! MAY NEED TO BE 
Changed for computers other than the coc 6000 series 


TOLC*(BOT*15.E-05)**2 
DO 6667 NUU=1,120 
Do 6667 NUT=1,120 
FvN(NUU,NUT)sO. 

6667 CONTINUE 

Do 308 NV=1.M 

CIR(NV.I)* 12.5663704 * ALP(NV) 

CIR(NV.2)= 12.5663704 

IF (PTEST.NE.O.) CIR(NV.2) = -1.0964155 * 0(NV) / BOT 
IF (QTEST.NE.O. ) CIR(Nv»2) a -1.0964155 * PV(NV) *BETA 

308 Continue 
I7Z=1 

NNVaTBLSCW(IZZ) 

DO 314 NV“1.M 
IZ=1 

NnN=T0LSCW(IZ) 

Do 316 NN=1.M 

APHI = ATAN(PHI (IZ> ) 

APSI = PSI (NN) 

XX=PV(NV)-PN(NN) $YY(1)=0(NV)-Q(NN) SYY (2) =0 (iNV) *0 (NN) 
ZZ=ZH(IZZ)-ZH(IZ) 

SnN = S(NN) 

Do 261 1=1.2 
YyY = YY(I) 

Call infsub (rot.fu(i) ,fv(D .fm<i) ) 

ApHI=-APHl tAPSl=-APSl 
261 CONTINUE 

IF (PTEST.NE.O. > GO TO 34? 

FvN(NV.NN)=FW(l)-FV (1 )*PhI ( IZ ) *FW (2) -FV (2) *PHI ( IZ) 

60 TO 312 

342 FvN(NV.NN)=FW(1)-FV(1)*PhI ( IZ ) -FW (2) *FV (2) *PHI ( IZ) 

312 IF (NN.LT.NNN .OR. NN.FO.M ) GO TO 316 
IZ=IZ*1 

NnN=NNN.TBLSCI»(IZ) 

316 Continue 

IF <NV.LT.NNV .OR. NV.EQ.M ) GO TO 314 
IZZ=IZZ*1 

NnV=NNV.TBLSCW(IZZ) 


M4TX 

IQ 

MATX 

2Q 

MATX 

30 

MATX 

40 

MATX 

50 

MATX 

60 

MATX 

70 

MATX 

80 

MATX 

90 

MATX 

100 

MATX 

llO 

MATX 

120 

MATX 

130 

MATX 

140 

MATX 

150 

MATX 

160 

MATX 

170 

MATX 

180 

MATX 

190 

MATX 

200 

MATX 

210 

MATX 

220 

MATX 

230 

MATX 

240 

MATX 

250 

MATX 

260 

MATX 

270 

MATX 

280 

MATX 

290 

MATX 

300 

MATX 

310 

MATX 

320 

MATX 

330 

MATX 

340 

MATX 

350 

MATX 

360 

MATX 

370 

MATX 

380 

MATX 

390 

MATX 

400 

MATX 

410 

MATX 

420 

MATX 

430 

MATX 

440 

MATX 

450 

MATX 

460 

MATX 

470 

MATX 

480 

MATX 

490 

MATX 

500 

MATX 

510 

MATX 

520 

MATX 

530 

MATX 

540 

MATX 

550 

MATX 

560 

MATX 

570 

MATX 

580 

MATX 

590 

MATX 

600 
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314 CONTINUE 

CALL MATlNV(FVN,M.CIRt2tDETERM.IPIV0T»INDEXtI20.ISCALE) 

IZZA = IZZ 
00 320 NZ=1»IZZA 
320 SECTRSTINZ) = 0. 

IZZ=1 

NnV>TBLSCW(IZZ> 

Do 614 NV*1*H 
IZ=1 

NNNsTBLSCH(IZ) 

VELIN * 0. 

Do 616 NN^l.H 

ApHI » ATANCPHKIZ)) 

APSI » PSI(NN) 

XX=PN(NV>-PN(NN) 

YY<H = 0(NV) - 0(NN) 

Yy(2) = 0<NV) ♦ Q<NN) 

ZZ“ZH(IZZ>-ZH(IZ) 

SNN * S(NN) 

Do 661 1=1.2 

YyY = YY(I) 

Call inesub (bot.fu(I) .fv(I) »fw(i) > 

aphi=-aphi 

APSI=-APSI 

661 Continue - (fv<1).fv(2)) • tak(aphi) >*cirinn.2) 

1 /FPI ♦ VELIN 

IF 1NN.LT.NNN .OR. NN.EQ.M > GO TO 616 
I7=1Z*1 

NnN=NNN*TBLSCWIIZ) 

616 CONTINUE 

CTCP = - (VELIN - 1, ) *2. * CIR(NV,2) 

SECTRST(IZZ) * SECTRSTd/Z) ♦ CTCP 
IF (NV.LT.NNV .OR. NV.EQ.M ) GO TO 614 
IZZ=IZZ*1 

NnV=NNV*TBLSCW(IZZ) 

614 CONTINUE 

return 

End 


HATX 610 
MATX 620 
HATX 630 
MATX 640 
MATX 650 
HATX 660 
HATX 670 
MATX 680 
MATX 690 
MATX 700 
MATX 71() 
MATX 720 
MATX 730 
MATX 740 
MATX 750 
MATX 760 
MATX 770 
MATX 780 
MATX 790 
MATX 800 
MATX 810 
MATX 820 
MATX 830 
MATX 840 
MATX 850 
MATX 860 
MATX 870 
MATX 880 
MATX 890 
MATX 900 
MATX 910 
MATX 920 
MATX 930 
MATX 940 
MATX 950 
MATX 960 
MATX 970 
MATX 980 
MATX 990 
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0VFRL4V(KlNr,Tt ,3.0) 

PPOGRAv aERODyM 

0 I mens t ON CPM(P) ,YCP(?) ,YY(?) .VOs'^l ,l!OU(l?0,?),FlH2),FV<?) * 

IXTLFGCFiO) .CHLFT (1?0 ,?) ,CLCC (1 ^ T| fG (FO) ,Fl DT (SO) ,CLA (2) ,SUM I 

?) ,AC(2) ,rH(2,S0) ,CCAV(?,50) .ClCL ( ?.f 0) , CP ( 1 20 ) , Fw < 2) 

3,DIFCIPS(?F) ,YIEGSV(?5) .ZIEGSV (2^) ,CLPT (120,2) , CL PR (120, 2) 

COMMON/ Al L/ P0T,M,RETA,PTFST,OTEST,THLSCW(S0),O(120),PN(12O)t 
1 PV (120) , ALP ( 120) ,S ( 1 PO) ,PSI ( 120) ,PHI ( 120) ,ZH (50) 

COMMON/TOTHRF/ C ip ( 120,2) ,SFCTPST (SO) 

COMMON/ONETHRF/TWIST (2) ,CREF , sPEF , CA VE,CLDES»STRUE» AR, ARTRUE, 

1 BTC0HT(2) .CONFIG, NSSWSV(2) ,Msv(2) ,KB0T, plan, IPLAN, MACH 

2 ,S5WWA(50) 

COMMON/THRECDI/SLOAD(3,50) 

COHHON/INSUB?3/APSI,APhI ,XX ,YYY,ZZ ,snn,tolcso 
1 FORMAT (/ 12X, *SECONO PLANFORM HORSESHOE VORTEX DESCRIPTIONS* / 

3 FORMAT (6F12. 5) 

4 FORMAT (IH1///58X,16HAER0DYNAMIC 0ATA///54X, •CONFIGURATION AFRO 170 

1NO.*F7.0 // ) AERO 180 

5 FORMAT (1H1,1RX*C0MPLETE C0NFIGURATI0N*31X*WING-B00Y CHARACTERISTICAERO 190 

IS* / 64X *LIFT* 9X *INDUCE0 DRAG (FAR FIELD SOLUTION)*// AERO 200 

2 16X A8 * CL COMPUTED ALPHA*19X *CL(WB)* 7X *CDI AT CL(WB)*AER0 210 

3 4X ,15HC0I/(CL(WB)**2) / 88X 1 2H ( 1/ (PI*AR) = F8.5 * )* ) AERO 220 

6 FORMAT (llX,2Fi5.5,15X,3F15.5) AERO 230 

7 F0RMAT(////4X,11H REF. CHORD, 6X,25HC AVERAGE TRUE AREA ,2X AERO 240 

1*REFERENCE AREA*9X*B/2* 8X,7HREF. AR,8X7HTRUE AR,4X , 1 IHmACH NUMPAERO 250 

AERO 260 
AERO 270 
AERO 280 
AERO 290 
AERO TOO 
AERO TIO 


AEPO 

10 

AERO 

20 

AERO 

30 

’AERO 

AO 

AERO 

50 

AERO 

60 

AERO 

70 

AERO 

80 

AERO 

90 

AERO 

100 

AERO 

110 

AERO 

120 

AERO 

130 

AERO 

lAO 

AERO 

150 

AERO 

160 


2ER/) 

8 FORMAT (8F15, 5) 

11 FORMAT (/// 47X ‘COMPLETE CONFIGURATION CHARACTERISTICS* // 

1 36X *CL ALPHA* 8X *CL (TWIST) ALPHA AT CL=0 Y CP Cm/CL 

2 CMO* / 27X *PER radian PER DEGREE* / 24X,7F12.5 ) 

12 FORMAT (//25X,*ADDITIONAL L0ADING*/24X*WITH CL BASED ON S(TRUE)* 

1 /67X34HL0AD DUE ADD. LOAD AT BASIC LOAD3X .27HSPAN LOAD AAERO 320 

2T SL COEF FR0M/8H STATION6X5H 2Y/B9X9H SL COEF ,4X8HCL RATIO, 4X7AER0 330 

3HC RATIO, 7X,14HT0 TWIST CL=,F9.5,3X,7HAT CL=05X,26HOESIREd CL AERO 340 
4 CHORD BD VOR/) AErO 350 

13 FORMAT (/ 47X, ‘CONTRIBUTION OF THE SECOND PLANFORM TO SPAN LOAD OAERO 360 


AERO 370 
AERO 380 
AERO 390 
AERO AOO 
AERO 410 
AERO 420 
AERO 430 


IISTRIBUTION* / ) 

15 F0RMAT(4X,I4,F12.5,5X,3F12.5,3X,3F12.5,3X,2F12.5) 

16 FORMAT (IHl) 

18 F0RMAT(////55X,21HTHIS CASE IS FINISHED) 

20 F0RMAT(///5X*DELTA CP TERMS FROM LE TIP TO TE TIP THEN INBOaRD 
IENDING with the TE OF ROOT CHORD *) 

21 FORMAT ( /54X*CMQ AND CLQ ARE COMPUTED*//) 

22 FORMAT(/38X*STATIC LONGITUDINAL AERODYNAMIC COEFFICIENTS ARE COMPijAERO 440 

ITED*//) afro 450 

23 FORMAT ( /59X*CLP IS COMPUTED*//) AERO 460 

24 F0RMAT(8F15.5) AERO 470 

25 FORMAT (/20X *X* IIX *X* IIX *Y* nx *Z* 12X *S* 5X *C/4 SWEEP* 4*AERO 480 

1 *DIHEORAL* 2X ‘LOCAL ALPHA* 2X ‘DELTA CP AT DESIRED* / AERO 490 

2 19X *C/4* 9X *3C/4* 42X *ANGlE*7X , ‘ANGLE* 4X,*IN RADIANS* 4X AERO SOO 

3 *CL =* FIO.S / ) AFRO SIO 

303 Format ( 1 2x ,9F 1 2 . 5 ) aero s20 

1013 F0RMAT(/47X*C0NTRIBUTI0N OF THE SECOND PLANFORM TO THE CHORD /iR OPAERO S30 
IAG FORCE*/) AERO S40 

1070 Format (//// sox, *induced drag, leading edge thrust and suctio’ aero s5o 

1 COEFFICIENT CHARACTERISTICS*/ AERO 560 

2 34X ‘COMPUTED AT ONE RADIAN ANGLE OF ATTACK FROM A NEAR FIELD SOl AERO 570 

3UTION* // aero 580 

4 S8X "SECTION COEFFICIENTS* 12X ‘CONTRIBUTIONS TO TOTAL CoEf.*/ AERO 590 

5 92X *FROM EACH SPANWISE ROW* / AERO 600 

6 38X *L. E. SWEEP* / AERO MO 
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7 15X 9X ^ ?Y/B* BX FX»C0T1 r./?H* SX 

ft sx c/^R« ftx*cnn« 9x lox ) 

1071 FORMAT (lOX ,T10^ BX, ftFl?,S) 

1072 FOPMAT (/// S7X,«TOTAL COFFF I C I FMTS ^ // ^ ^ 

1 36X IPHCOII/Cl = FlO.5 fSX = * 

1074 FORMAT(////10X<^1NVALID LEADING EDGE ’ 

1* G14.S * FOP THE* IS * SPANWiSE POW. CS IS WRONG. ) 

4445 F0RMAT(//////////S6X,4HCLP=»F9.5////) c/y//\ 

4446 FORMAT (//////////42X»4HCMQ=»FR*5, 1 OX *4HCLQ=*FS.5////) 


THF COST^^P 


C 

c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 


PART 3 - COMPUTE OUTPUT TERMS 


RAD = S7.2957G 

TWST ’= TWIST(I) ♦ TWIST(2> 
ALREF = I 


THE TOLERANCE SET AT THIS POINT IN THE PROGRAM M*Y NEED TO RE 
CHANGER FOR COMPUTERS OTHER THAN THE COC 6000 SERIES 


TOLC= *0100*807 
TOLC50 = TOLC*TOLC 
QINF=1 . 

NSSW^NSSWSV ( 1 ) ♦NSSW5V(2) 

IE (RTCDHT ( 1 ) *NE*RTCDHT (2) > GO TO 794 
SUMPHI=0 
DO 801 J=1»NSSW 
801 SUMPHI=SUMPHI^ABS(PHI (J) > 

IF (SUMPHI.EQ.O.) GO TO 921 


part 3 - SECTION 1 

COMPUTE lift and PITCHING MOMENT FOR WINGS WITH DIHEDRAL 


GEOMETRY FOR TIP TRAILING LEGS 


794 CPM(1)=CPM(?)=YCP(1)=YCP(2)=IM=CLT=CLNT=NSSWI=0 

NSSW2 = NSSW3 = NSSWSV(l) 5L=I 

NSCW = MSV(l) / NSSWSVU) 

60 TO 798 

796 NSSWI = NSSWSV(I) i 

NSSW2 » NSSW SNSSW3=NSSWSV(2)$L=NSSWSV(1»*I 

NSCW = MSV(?) / NSSWSV(?) 

798 I = IM ♦ I 

J = IM ♦ 2 


iuy=2 


DIEFCRI=0. 

APHI=ATAN(PHI (I) ) 
TLX1»PN(I)-5(I)*TAN(PSI (I) ) 
TLX2=PH( J) -S( J) *TAN(PSI ( J) ) 


CLFTL6=TLX1-TLX2 
XTLEOn >=TLXl/2.*TLX2/2. 

YLE6=0 ( I ) -S ( I ) *COS ( APH I ) 

IE (NSSWI .EO.O) YLEGSVI 1)=YLEG 
ZLEG*ZH ( I ) -S (I ) *SIN ( APHI ) 

lE(NSSWl.EQ.O) ZLE6SV( 1 ) =ZLEG 


AERO •' ?n 

Arps ' 10 

APPO ^40 
AFRO 

AFRO 860 

-AFRO ‘^Gl 
AFRO ^6? 
AFRO ^70 

afro ^80 

afro 890 
AERO 700 
AERO 710 
AERO 720 
AFRO 730 
AERO 740 
AERO 750 
AERO 760 
AERO 770 
AERO 780 
AERO 790 
AERO ftOO 
AFRO ftlO 
AERO "^20 
AERO ft30 
AERO 840 
AERO 850 
AFRO 860 
AERO 870 
AERO 880 
AERO «90 
AERO 900 
AERO 910 
AERO 920 
AERO 930 
AERO 940 
aero 950 
AERO 960 
AERO 970 
AFRO 980 
AFRO 990 
AEROIOOO 
AEROIOIO 
AER01020 
AEPO1030 
AER01040 
AER01050 
AER01060 
AER01070 
AER01080 
AER01090 
AEROllOO 
AEROniO 
AERO1120 
AERO1130 
AEROn40 
AER01150 
AEROl 160 
AEROlWO 
AERO1180 
AER01190 
AER01200 
AEROIOIO 
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rF(NSSWl .FO.NSSWSVf U ) GO TO 8SC 
Gf) TO 857 

rt5o nn SOSO tt=i*l 

TF( { ARS ( VI FGSV( m-YL EG) . LT.TCLC ) .AND. C ARS( 7L EGSVC in-ZLEG).LT, 
ID) iirFFCRl=niFCrRS( IT) 

505D rONTTNUF 
857 nn 807 KV=7*NSCW 
NVT=NV-1 

807 XTI FGI NV)-XTJ.FG(NVT)“CLFTLG 
NCTI =0 tNA^l $NB=NSCW 
80 T nn 87i NV=NA*NR 

VniH NV . 1 ) = Vnn|NV.7 ) = UCU(NV , l )=UCU(NV.7)=0. 

nn 809 NN=|,M 

I7^(NN-1 )/N5CW+l 

APMI^ATANCPNI (I 7 m 

AP5T = PS H NN) 

XX^XTJ FGINV)-PN(NN) 

Yvn ) = Y| FG-OC NN) 

YYI7I=YI FGfO(NN) 

77 = 71 FG 

5NN = 5(NN) 

nn 877 1=1.7 

YYY = YYCI) 

CAM fNFSUB IHOT.FliC I ) ,FV( n ,Frt( I) ) 

APHr=-APHI SAPS 1 = -APS I 

877 CONTINUF 


9 nni on 809 Txx=i .7 

nniH Nv/. I XX )=unU(NV t 1 XX) + ( ( FU ( l ) ♦Fiji 7 ) ) ♦CI R ( nn. I XX ) ) / l 7. 56637 l 
VntJI NV. IXX J =vnu CNV.tXX) + |(FV(l)+FV(7))*CIR<NN.IXX)l/l7. 566 3 7 1 
CONT r NUF 
NCTI =NCTI 

IF INCTI-7) 810.811,812 


809 
87 3 


C 

C 

C 


GFOMFTRY FOR SPANWTSE BOUND VORTICES 

8in NA=NSCW+1 
N8=7*NSCW 
JA= T M*NSCW4-1 
Yl FG=0C JA) 

71 FG = 7HI TM4-1 | 

nn 818 .i=i.Nscw 

.IK=r P + NSCW + J 

nv=.unscw 

818 XTI FG( NV) = PN(JK ) 

r,n TC 80 3 

GFDMFTRY Al DNG RIGHT TRAILING LEGS 

811 NA=7«NSCVi + l 
NR=3*NSCW 
niFFCR7=0. 

.IK= I 84tNSCW+l 

APHI =ATAN( PHI ( I M+1 ) ) 

Yl FG = 0(JKM-S( JK)»COS| APHI ) 

IFINSS-Jl .F0,0) YLFGSVI IUU)=YLFS 
71 FG=7H« IM<^1)^S(JK)*SINIAPHI) 

IF (NSS^l .FC.O) 71 FGSV( I bU) =ZLEG 
Tl XI =PN(JKl+S(JKlfrTAN(PSI( JK) ) 

JK=JK+1 


AERni220 
AERC1230 
AER01240 
TCLAERC1250 
AERG1260 
AER01270 
AER01280 
AER01290 
AFR01300 
AER01310 
AER01320 
AER01330 
AER013A0 
AER01350 
AERC1363 
A ERG1370 
AER01380 
AER01390 
AtKQlAOO 
AEROIAIO 
AERC1420 
AER01430 
AER014A0 
AER01450 
AER01A60 
AER01470 
AER01A80 
AERCIA90 
AER01500 
AERC1510 
AERQ1520 
AER01530 
AER015A0 
AER01550 
AER01560 
AFRC1570 
AER01580 
AFR01590 
AER01600 
AERG1610 
AERC1620 
AE301630 
AER016A0 
AERG1650 
AtR0I660 
AER01670 
AERD1680 
AERQ1690 
AERG1700 
A6Rni7lO 
AER01720 
AERG1730 
AER01740 
AERC1750 
A6R0I760 
AERG1770 
AERC1780 
AER01790 
AERG1800 
AFR01810 
AER01820 
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TLX7=PN(JKI^SIJK)*TAN(PSI( JK)) 

CRTTLG=TlXl-TtX? 

XTLFGfNAI=TLXl/2.<-TLX2/2. 

NAA^NA^-l 

IFINSSWl.FO.NSSWSVC n I GO TO 851 
GO TO 853 

IF* ?ABSC YLEGSVf ITI-YLEG) . LT.TOLC ) . AND. ( AftS< 7L EGSV U T » -ZLEG ) 
ICI) DIFFCR2=DIFCIRS( IT) 

5051 CONTINUF 
853 DO 810 NV=NAA»NB 
NVT=NV-I 

810 XT! FG(NV)-XTIFG( NVTI-CRTTLG 
GO TO 803 

COMPUTE LIFT AND PITCHING MOMENT FOR FACH ELEMENTAL PANEL 


C 

c 

c 


835 


834 


B17 YYI 1 )=YYf 21=0 

IF < TM-NF.NSSWl ) GO TO 834 
on 835 IXX=U2 
DIFCIR=DIFFCRI 
DO 835 NPOS=l»NSCW 
DIFCIR*OlFCIR^CIR<NPOSf IXX) 

CON=l. 

IF (NPOS.EO.NSCW) C0N=.75 

CHI FT(NPnS*IXXl=CLFTLG*CCN*OIFCIR4VOU(NPOS, IXX)*(2./SREF) 

CLPTC NP05*IXX) = CHLFT< NPOS, IXX ) ♦ I 0 ( NPOS ) -S I NPOS > ) *2. 

CONTINUF 

IFIN5SW1 *F0«O) DIFCIRSI 1 )=OIFCIR 
DO 815 IXX=l»? 

D!FCIR=DIEFCR? 

DO 815 NPOS=l.NSCW 
jKsIM^NSCW+NPOS 
Jl =1 IM+1)4NSCW^NP0S 
JMsNSCM^NPOS 

jn=?*nscw^npqs 

IF I IH,F0.fNSSW2-l ) ) GO TO 836 

DIFCIR=DIFCIR>CIRC JL,IXX)-CIft(JK,IXX) 

CON=l. 

IF { NPnS.FO^NSCW) CON=.75 

CHIFTI JL,IXXI=CRTTLG*CON»OIFCIR*VOU( JN, IXX)^(2./SRFF ) 

Cl CC UK* 1XX) = < 2. /SRFF) 40181 JK, I XX) *2. ♦SIJK)*Cn<;(APHn* < I 

UXXI+VOUUMf IXX )4TANIPSI(JK) ) ) 
CLPBUK*IXXI=CLCC(JK,IXX)*Q(JK)*2. 
CLPTUL*IXX)=CHLFTCJL,IXX)4IQUK)+SUK) )42. 

YYI lXXI=YYf IXXI4CCLCCIJK. I XX ) 4C HLFT UK • IXX) 1*2. 

CPMI IXXI=CPMI IXXl+ICLCCUKf IXX)4XTLEGUM)48PTA>CHl FTUK*IXX 
1 C NPnS)*RFTA»*2./CREF 

YCPI IXX)=YCPf IXX)<“fCLCCUK,IXXl*QUK)4CHLFTUK» IXX) 4I 0(JK)- 
icnsf APHini/BOT 
815 CONTINUE 

IFINSSWI.FO^O) DIFCIRSI IUU)=DIFCIR 
CLT=CLT+YY( 1) 

CLNT»CLNT>YYI2) 

IM=IH-H 

IF! N5SUI.E0.0) IUU=IM'i-2 
IFIIM.FO.NSSWSVm) CLWNGT=CLT 
IF! IM.FO.NSSWSVI 1) ) CLWING=CLNT 
IF I IM.GF.NSSW2I GO TO 816 

NCTI =*1 

DO 817 IXX=l*2 


836 


AERG183C 
AER01840 
AER01850 
A6RC1860 
AERQ1870 
AER01880 
AER01890 
,LT.TCLAER01900 
AER01910 
AER01920 
AER01930 
AER01940 
AERC1950 
AER01960 
AERC1970 
AER01980 
AER01990 
AER02000 
AER02010 
AER02020 
d ERC2030 
AERQ2040 
AEP02050 
AER02060 
AER02070 
AER02080 
AER02090 
AER02100 
AER02110 
AERQ2120 
AER02130 
AERC2140 
AER02I50 
AER02160 
AERQ2170 
AER02180 
AERG2190 
AER02200 
AER02210 
AER02220 
AER02230 
UCU( AER022^0 
AER02250 
AER02260 
AER02270 
AER02280 
)4XTLEGAERC2290 
AER 02300 
SUKI4 AER02310 
AER02320 
AER02330 
AERO23A0 
AER02350 
AER02360 
AERC2370 
AERC2380 
AERQ2390 
AER02400 
AER02410 
AER02A20 
AERG2430 
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on «17 NV -l.NSCW 
MY=NV4'P*NSCW 
XT! FOCNVl-XTI FG(NYJ 
fll7 VniMNV. IXX»=VnU(NY. TXXJ 
GO TO 810 
C 

r. SUM I TFT AND PITCHING MCMENT FOR ENTIRE WING 

C 

81 A YY( 1 ) = CI T*SRFF/STKUF 
YY(7)=n NT^SRFF/STRUE 
NIJP = NSSW3 ♦ 1 
YTI FG(NUPI=0. 

XTI FG{NUP)-0 
TN0=1 

IF ITWST .FO.O.I IN0^2 
on 817 TXX=IN0*7 
on 870 JS8W=L.NSSW7 
OAO( I XX, JSSWI=0 
SI 0T( JSSW) = 0 
APHr^ATANCPHICJSSWI 1 
J1 =( JSSW-l l*NSCW+l 
K=JSSW-I+1 

870 YTIFG( K MO lJU - S ( JL) ♦CC S ( APH I) 
on 817 TNC = 1,NSCW 
on 818 JNS=L,NSSW? 

JK=( JNS-l MNSCW+INC 
K = JNS~I ♦! 

818 XTlFGf K MCHLFTCJK, IXX) 
on 817 INS=L,NSSW7 
JK=f INS-1 MNSCW+INC 
APHI^ATANIPHT C INSI ) 

CAM FTIUP <OUK),CHTLF,*UNUP,YTLFG,XTI FGI 
T= SRFF/(7.4‘SC JKMCnSI APHI ) ♦CAVE I 
Si OTI INS) = SLDT( INS)*CHTLF^T 

CICCIJK.IXX) = ICLCC( JK, IXX) ♦ CHTl F I ♦ T 
817 SI OADI I XX, TNSMSLnADC IXX, INS1+ CLCCUK,IXX| 

IF (IN,NF,NSSW) GO TO 796 

Cl AI7MCI NT /At REF 

CMCI =CPMC7)/CLNT 

CMn=CPMM )-CMCL^Cl T 

YCP(7MYCPI7)/(CLNT/2,) 

on 840 T=1,NSSW 

SI OTIIMSI DTI I)/YYC?1 

IF ITWST •FO.0,1 SLOA0(l,n=0. 

IF ITWST •NF.O,> SLOAOI I, IMSLOAOI I, I )/YY( 11 
840 SI 040(7,1) = SI nA0(2, II/YY(2) 

CRI =0* 

on 860 IAM=l,M 

860 CRI =CRI +CLP8I IAM,7)+CLPT( IAN, 2) 

Cl P=CR| /C,08778+7.*B0T) 

GO TO 901 
C 

C PART 3 - SECTION 7 

C COMPUTF I IFT ANO PITCHING MCMENT FOR WINGS WITHOUT DIHEDRAL 

C 

971 on 901 NV=1,7 
SUM! NVI=0 
on 901 I=l,M 

SilM(NVI=SUMINV)+CIR( 1 ,NV)^S( I) 

IF INV.FO,UAND,I.EO,MSV(1 ) 1 CLWNGT = SUMIIMB. / SREF 

IF INV.F0.7,AN0.I,E0.MSV(1 ) ) CLWING = SUMI2)^a. / SREF 


AER02440 

AER02450 

AER02460 

AER02470 

AERQ2480 

AERC2490 

AER02500 

AER0251O 

AER02520 

AER02530 

AER 02540 

AER02550 

AER02560 

AER02570 

AER02580 

AER02590 

AER02600 

AER02610 

AER02620 

AER02630 

AER02640 

AERC2650 

AER02660 

AER02670 

AER02680 

AERG2690 

AER02700 

AERG2710 

AER02720 

AERG2730 

A6R02740 

AER02750 

AER02760 

AERQ2770 

AER02780 

AER02790 

AER02800 

AERG2810 

AEPG2820 

AERG2830 

AERG2840 

AERG2850 

AERG2860 

AER02870 

AER02880 

AERG2890 

AER02900 

AER02910 

AER02920 

AER02930 

AER02940 

AERG2950 

AERG2960 

AERQ2970 

AERO290O 

AER02990 

AERG3000 

AERQ3010 

AER03020 

AER03030 

AER03040 
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901 CONTINIIF 

Cl T = SUM(l>/SRFF 

CLMT = SUMI7I/SREF 

IF (KBOT.FO.ll 
CI-WNOT = CLT - CLWNOT 
Cl WING = CLNT- CLWING 
BOO CRI = 0, 

on 90S T = 1 

CRI =CR« +( 0 I n^CIRl I *2 )*2.*S( n I ♦?. 
C* CCI 1*1 1=C IRC I .11*2. /CAVE 
90S ClfXII*7>=CIRCl*2)*2./CAVe 

COMPUTF CLP 


AERC3050 
AERC3060 
6FRD3070 
AERU3080 
AERQ3090 
AERG3100 
AER031 10 
4ERL33120 
AERQ3130 
AERC31^0 
AERa3150 
AFRG3160 
AErt03170 
AFR03130 


C 

c 

c 

c 


Cl P= r.Rl / C SR FF*Bni ♦0.00725) 

Cl AC 21=CLNT 
on 922 IXX.= l*2 
SA=SR=Sr=0* 

I ^ 0 

on 920 JSSW=1*NSSW 

SI nTCJssw)=o 

SI OAOI !XX* JSSWI=:0 

NSCW = TBLSCWCJSSW) 

on 9?0 JSCW=l.NSCW 

IFITWST •F0*0. . AND. IXX.EQ. 1 ) GO TO 930 
I =1^-1 

SA==SA+C1R1 I.TXXl^SC I ) 

SB=sa*c IRC I. ixxi^oc n*sc n 

SC=SC^C IRC I *IXX )^PN( I )*SII )*13ETA 

SI OAOC 1 XX* JSSWI = St OADl IXXf JSSW) ♦ I anT*CT R ( I 

GO TO 920 

930 SLOADIl *JSSWI=0. 

920 CONTINUF 

IFITWST •F0*0. • ANO . 1 XX. EO . 1 > GO TO 932 
yCPC IXX)=SB/CSA^ROTI 
ACC IXX)=SC/CSA^CRFF) 

GO Tn 922 

932 YCPCl 1-ACC11=0. 

922 CONTINUE 

CMCI^ACI2) 

CMO=l ACC I l-AC( 2))*CLT 


IXX) )/(7.^SUrtC IXX) ) 


PART 3 - SECnCN 3 

COMPUTE AND PRINT FINAL OUTPUT DATA FOR ALL WINGS 


AERC3190 

AtR03200 

AER03210 

AERG2220 

AtRO'3230 

AEHC3240 

AERa32SQ 

AERC3260 

AER03270 

AERU3280 

AERC3290 

AER03300 

AER03310 

AERC3320 

A6RQ3330 

A ER03340 

AERQ3350 

AERC3360 

AER03370 

AER03380 

AER03390 

AERC3400 

AER03410 

AERC3420 

AERC3430 

AER03440 

AERC3450 

AEBG3460 

A ERG3470 

AER03480 

AtRC3490 


903 no 902 IXX=1.2 

JN = 0 

DO 902 JSSW=1*NSSW 
CH (IXX.JSSWI=0 
NSCW = TBLSCWCJSSW) 

DO 904 JSCW=UNSCW 

CH I I XX* J SSW 1 = ( -2 • 0 ) ♦ C P V I JN)- PN C JN ) ) *BF T A+CH (IXX.JSSW) 

904 CONTINUF 

CCAVC IXX* JSSWI =CHC IXX»JSSW)/CAVE 

Cl Cl C I XX* JSSWI = Sin ADC IXXf JSSWl/CCAV C I XX* JSSW) 

902 CONTINUF 
CLD=CLDFS 

IFCCI DFS-FO-ll ) Cl 0=1* 

CPCn^^ I CCrCCCI.n^CLCCC If2)*(CLD -CLT)/CLNT)*CAVF/C?.*(PN( I )- 


AER03500 

AERC3510 

AFR03520 

AFRC3530 

AER03S40 

AER03550 

AFRG3560 

AERG3S70 

AER03580 

AERC3590 

AER03600 

AtR036t0 

AER03620 

AERC3630 

AER03640 

AEPG365C 
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1 PV(H ) * BETA > 

lOZO CONTINUE 

WRITE (6,4) CONFIG 
If ( PTEST.NE.O. ) 

IF ( QTEST.NE.O. 


) 


WRITE (6,23) 
WRITE (6,21) 
WRITE (6,22) 


head * OH 


IF ( PTEST.EO.O. .AND. OTEST.EQ.O, ) 

Write (6,25) cld 

head » 8HDESIRED 
IF (CLOES.EQ.lI. ) 
lEND » 11 

IFICLOES.NE.II.) IENO=l 
Do 5000 IUTK=1,IEND 

IF(IEND.EO.II) CLDES»(FLOAT(IUTK)-1.)/10. 
IFICLOES.EO.O.) CL0ES=-.) 

Nr a 0 

DO 3006 NV»l,NSSW 
N<;CW » TBLSCW(NV) 

Np X nr « 1 

nr » nr , NSCW 

PHIPR X ATAN (PHI (NV) ) * 

Si,OADn,NV)=0. 

IF (NV.EO. (NSSWSV(l) 

Do 3006 ixNP,NR 
IF ( lUTK.GT.l ) 

PnPR X PN(I) • beta 
PVPR = PV(I) * BETA 
PSIPR = PSKl)* RAD 

(6,303) PNPR,PVPP,Q(I) ,ZH(mv) ,S(I) ,PSIPR, PHIPR, 
3006 SlOAD(3,NV) =SlOAD( 3,NV) ♦0LCC(I,2)*r(.DES/CLNT*CLCC(I,l 
ILT/CLNT 

IF(IUTK.gT.I) go to 3007 
Write (6,7) 

■»««-, crff,cave,strue,sref, BCT,AR,APTRLE,MACH 

Jv07 CONTINUt 


IP(PTE<;T.N€:.0.>WRITF(6,4445) CLP 
IF(PTEST.NE*0*) GO TO 44^4 


RAO 


1> > WRITE I6*l) 


GO TO 3006 


AER03660 
AERO3670 
AER03680 
AER03690 
AER03700 
AER03710 
AERO3720 
AEPO3730 
AER037A0 
AER03750 
AER03760 
AERO3770 
AER03780 
AEPC3790 
AERC3800 
AERO301O 
AER03820 
AERC3830 
AEP03840 
AER03850 
AER03fl60 
AERO307O 
AER03880 
AER03890 
AEPC3R00 
AER03910 
AER03920 
ALP(I),CP(I) AER03930 
) -CLCC a *C AEPC3940 
AERC3950 
AER03960 
AER03970 
AER03980 
AEPC3990 
AERC4000 
AER04010 
AERC4020 
AEPC4030 


COMPUTE CMOtCLO 

CmQ» 2,0*CMCL*CLNT/(0,0R7?S*CPEP') 
Ci,Q»2.O«CLNT/(O.O072S*CRrF> 

IF (OTEST.NE.O. ) WR t TF (6ta446 ) CMQ^CIO 
IF<OTEST.NE,0.) GO TO A4a4 

COMPUTE INDUCED ORoG 

NSV»NSSWSV(U ♦! 

HtOT=HSV(1>*1 

If (KBOT.FO. 1 ) (50 TO 1001 

NSV»NSV^N5SWSV (?) 

MtOT=MTOT^MSV(?) 

1001 Call COTCLS (AR*ARTRUFfNSSWSV<K'0OT) ,MTOTfNcVfCDI*CDIT) 
ClAPOsClA (?)/ 87. 29578 
AjLPOa- (CLT/Cla ( 2H *57, 29570 

alpd=cldfs/clapo*alpo 

AlPW=1,/CLAPO 

ClW8=CLWING*ALPD/57.2957r^CLWNGT 
CpiWS = CDI /(CLWfl*CLWR) 

IF (lUTK.EO.n WRITE (6,5) HFAO,CDIT 


AERC4040 

AERC4n50 

AERC4n60 

AFR04070 

AERC4O0O 

AFRC4090 

AFP(!4100 

AFR04110 

AFRC4120 

AFRC4130 

AFR04140 

AERC4I50 

AEPC4160 

AFRC4170 

AERC4180 

AER(}4190 

AER04200 

AFR04210 

AFP04220 

AEPC4230 

AEP04240 

AFRC4P50 

AEPC4260 
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APpf)4?70 

5000 WRITE (6,6) CLDES. ALPD.CLWB.COI.COIWB , AEROA?80 

WRITE(6,ll) CLA(2).CLAPD.CLT,ALPO,YCP(2).CMCL,CMO AER04290 

WRlTE(6fl2) CLT AER04100 

NR « J = 0 AER04110 

DO 1004 NV*lfNSSW AER04'^20 

BCLCC»BADLAE*BASLO*0* AERO4330 

NSCW = T0LSCW(NVI AER04340 

NP s NR ♦ 1 AERO4350 

NR s NR ♦ NSCW AERO4360 

DO 1002 I»NPfNR AER04370 

ADLAE*CLCC ( I *2) *CLT/CLNT AERO430O 

BSLDaCLCC < I * 1 ) -AOt AE AERO4390 

BCLCC=BCLCC^ClCC (1*1) AER04400 

BADLAEsBADLAE^ADLAE AER04410 

BASLDaBAStD^BSLD AERO4420 

1002 CONTINUE AERO4430 

j * J ♦ NSCW AER04440 

YQ = Q(J) / BOT uofTPiiL AERO4450 

I BASLDtSL0ADI3.NV),SL0T<NV) AEROA40O 

WRITE (6tI070) AEROAA90 

CTHRUST a CSUCT a CDRAG »0. AER04S00 

NN*1 AEROA510 

DO 1050 NV=1»NSSW «„-t» AERn4S20 

IlSSio : • CL..3. / .STPUE • 4. ■ BC„ ..Pn4P3. 


- SSCTRST 

r COS ( AT AN (SSWWA(NV))) 
r SSCTRST / CSSWWA 

IF (NV.EO.l) GO TO 1060 
NN a NN ♦ TBLSCW(NV-l) 


CSSWWA 

SSCSUCT 


1060 PHIPR 
CORAGS 
CDRAG 
CTHRUS5 
CTHRUST 
CSUCTS 


ATAN IPHI(NV)) 

SSC0RAG*A. *B0T*2.*S(NN)*C0S (PHIPR) /SREF 

CDRAG ♦ 2»0 * CORAGS 

SECTRST(NV)*2.*S(NN)*GOS(PHIPR) / SRFF 
CTHRUST ♦ 2.0 * CTHRUSS 

IF^THE ABSOLUTE^VALUE^Of'^THE LEADING EDGE SWEEP ANGLE IS GREATER 

J^aJ 80 ?IgreJs no suc?ion contribution is computed 
ir f r^^wWA .1 0#17365 ) CSUCTS * 0* 

IF ( CSSWWA .LT. 0. ) WRITE <6. IOTA) CSSWWA, NV 

CSUCT a CSUCT ♦ 2.0 * CSUCTS 

swale a ATAN(SSWWA(NV)> * RAD 

YQ » Q(NN)/ BOT 

IF(NV.E0.<NSSWSV(1)*1)> write (6, 10 13) 


AEROASAO 

AEROA550 

AEROA560 

AER0A570 

AEROAS80 

AEROA590 

AEROA600 

AEROA610 

AER0A620 

AEROA630 

AER0A6A0 

AER0A6A1 

AER0A6A2 

AER0A6A3 

AER0A6AA 

AEROA650 

AEROA660 

AEROA670 

AEROA680 


1050 ij;?^r6:iS7?rN:irolUArE:rsjD;iG.SSCTRST,SSCSUCT.CDRAGS,CTHRUSS^ 

^ r'ltr'T e- 


^ CSUCTS 

CORA6P » CDRAG / (CLA (2> *CLA (2> ) 
write (6fl072) CORAGPfCTHRUSTfCSUCT 
4444 WRITE<6*18) 

WRITE(6.16) 

RETURN 

end 


AER04700 

AER04710 

AER04720 

AER04730 

AER04740 

AERO4750 

AER04760 
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5URR0UTTME CnrrL«^ ( AR»/^RTPUe: • TsEMSPf MTOT»NSV,CDI ,CniT) 
dimension ETAN (51 ) ,GA'^PR (51 » 1 ) *FTA (41) * GAMMA (41 ) *VF (41 ) fP (4] ) , 
lFVN(41t41) 

COMMON/ all/ RnT*M,BETA,PTFST,OTrSTtTRLSCW (SO) »O(i;?0)tPN(l?0), 

1 PV (1?0> » ALP (IPO) tS ( IPO) ,PSI ( IPO) »PHI ( 1?0) t/HlSO) 

COMMON/THPECOI/SLOAD n,50) 
no 15 T=lt41 
DO 15 .1=1 f41 
15 FVN(I*J)=0 
SPAN=P.*ROT 
CAV0 =SPAN/ARTpUE 
PI=.314159?65F^0l 
NST=I5FMSP^l 
NN=MTOT 

DO 101 N=1*ISFMSP 
NM=NSV - N 
NSCW=TPLSCW(K)M) 

NN=NN-NSCW 

ETAN(N) =ASIN(-0(NN) <^?./SPAN) 

GAMPR(Mf I ) =SLOAD ( 3 , NM ) <^CAVB / { ? . ^*SPA^ ) 

101 CONTiNtlE 
ETAN(NST>= PI/?. 

GAMPR(NST» 1 ) =0 
DO 7 NP= 1*41 
ANP=NP 

7 ETA(NP)= ( ANP-Pl . ) »Pl/4?, 

DO 102 JK=?lf4l 

CALL FTLUP (ETA ( JK) ^CtAMMA ( JK) , 1 t MST,ETAN»GAMPR) 

102 CONTINUE 

DO 600 NY=P?#41 
ETA(NY)=SIN(ETA(NY) ) 

NR=4?-Ny 

ETA (NR) =-ETA (NY) 

600 GAMMA (NR) =GAmmA (NY) 

DO 5B9 NU=21*41 
ANU=NU 

DO 14 M=1 ,41 
AN=N 

NNUD=I A0S (N-NIJ) 

VE(N)=COS( ( (AN-?l.)*Pr>/4?,) 

IE (NNUn.NE.O) GO TO 9 

B(N)=(42.)/(4.0^COS( ( ( ANU-2 1 . ) «P T) /4? . ) ) 

GO TO 14 

9 IF(MOD(NNUD,?> .EQ.O) GO TO 12 

B(N)=VF(M)/( (42,)*(FTA(N)-ETA(NU) )»*2) 

GO TO 14 
12 B(N)=0.0 
14 CONTINUE 

DO 5B9 NP=?1,41 
NUST =TARS(NU-P1) 

IF(NUST,FQ*0) GO TO 589 

IE (M0D(NU5T,2> .EQ.O) GO TO 589 

NPST=IABS(NP-PO) 

IF(MOD(NPST,2) .EO.O) GO TO 589 
NPNU0=1ARS(NP-NU) 

IF(NPNUD,EQ.0) GO TO 589 
IF (MOD(NPNUO,?) .EQ.O) GO TO 589 
FVN(NU*NP)=2.0<»B{NP)/?l.*COS( ( ANU-21 . ) "PI /4? . ) 

IT=42-NU 


CDIC 

10 

CDTC 

20 

CDIC 

30 

CDIC 

40 

CDIC 

50 

CDIC 

60 

CDIC 

70 

CDIC 

80 

CDIC 

90 

CDTC 

1 00 

CDIC 

110 

CDTC 

120 

CDIC 

130 

CDIC 

140 

CDIC 

150 

CDIC 

160 

CDTC 

170 

CDTC 

1 80 

CDIC 

190 

CDIC 

?00 

CDIC 

Pio 

r.DTC 

P20 

CDIC 

^30 

CDTC 

240 

CDIC 

^50 

CDTC 

26 0 

CDIC 

270 

CDIC 

280 

COIC 

^90 

CDIC 

^00 

CDTC 

MO 

CDIC 

120 

COIC 

130 

CDTC 

140 

CDIC 

150 

COIC 

160 

CDIC 

170 

COIC 

180 

CDIC 

190 

COIC 

400 

COIC 

410 

CDIC 

^20 

CDIC 

410 

COIC 

440 

CDIC 

450 

CDIC 

'^60 

CDIC 

-.70 

CDIC 

480 

CDTC 

490 

CDIC 

ROO 

CDTC 

^10 

COIC 

«=^20 

CDIC 

s30 

COIC 

S40 

CDTC 

550 

COIC 

560 

CDIC 

570 

CDIC 

580 

CDTC 

590 

CDIC 

500 
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FVN(NU*ITT)=?*0«R(TTT) /21*^^C0Sl 
FVN ( IT f NP) =FVN (MU» T TT) 

FVM (IT.TTT) =FWN(NU»MP) 

SH9 CONTINUE 

CCC = 0.^'- 
DO 10 N=lt41 

I 0 CCC = CCC^ (OA'^MA (N) ^GAMMA (N) ) 

CCD=0.0 

00 11 ''TJP^l.41 
DO 11 ^!=1»A1 

CCn=CCD-?.0»FVN(NUP^N) (DAMMA (NUP) *GAMMA (N) ) 

II continue 

cni=Pi<^AP/4.« (ccc^cco) 

COIT=l -/(PT^AP) 

RETURN 

END 


CDIC 

CDTC 

COTC 

COTC 

CDIC 

CDIC 

CDIC 

CDTC 

CDIC 

CDTC 

CDIC 

CDIC 

CDIC 

CDIC 

CDIC 

CDIC 

COIC 

CDIC 


in 

" ?0 
^30 

^50 

<S60 

<s 70 

*^90 

700 

710 

7?0 

730 

7A0 

7S0 

76 H 

770 

7R0 
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SUBROUTINE f«ATINV< A t , M,DETERM,I PI VOT, INDEX ,NHAX, I SC ALE ) 

DOCUMENT DATE 08-01-68 SUBROUTINE REVISED 08-01-68 

MATRIX INVERSION WITH ACCOMPANYING SOLUTION OF LINEAR EQUATIONS 

DIMENSION IPIVOT(N» ,A(NMAX»N) rB<NMAX,H> , INDEX(NMAX,2 ) 
EQUIVALENCE (IRCW,JROW|, ( ICOLUM , JCOLUM I , (AMAX» T, SWAP) 

INITIALIZATION 

5 1SCALE*0 

6 9-1 = 10.0**100 

7 R2=i.0/Rl 
10 0ETERM=1.0 
15 DO 20 J=1,N 
20 IPIV0T(J)=0 
33 DO 550 1=1, N 

SEARCH FOR PIVOT ELEMENT 

40 AMAX=0.0 


45 

DO 

105 

J=1,N 




50 

IF 

( IPI 

V0T( J)- 

-11 

60, 

105, 60 

60 

DO 

100 

K = 1,N 




70 

IF 

( IPI 

VOTCKJ- 

-1) 

80, 

100, 740 

83 

IF 

( ABS 

C AMAX)- 

-ABSIAI J 

,K ) 1)85, 


85 1R0W=J 
93 ICOLUM=K 
95 AMAX=A(J,K) 

103 CONTINUE 

105 CONTINUE 

IF (AMAX) 110,106,110 

106 0ETERM=0.0 
ISCALE*0 
GO TO 740 

113 IPIVOT( !COLUM)=IPIVOT( IC0LUM)+1 

INTERCHANGE ROWS TO PUT PIVOT ELEMENT ON DIAGONAL 

130 IF ( IROW-ICOLUM) 140, 260, 140 

140 DETERM=-DETERM 

150 DO 200 L=1,N 

160 SWAP=A( IROW,U 

173 A(IROW,L)=A< ICGLUM,L) 

200 A( ICOLUM,U=SWAP 
205 IF<M) 260, 260, 210 
210 DO 250 L*l, M 
223 SWAP*B( IROW,L) 

230 IROW,L)=B< ICCLUM,L) 

253 8( ICOLUM,L)*SWAP 
263 INDEXC 1,1 J=IPOW 
270 INOEX< I ,2)=IC0LUM 
313 PIVOT=A(ICOLUM, ICOLUM) 

IF IPIVOT) loco, 106, 1000 

SCALE THE DETERMINANT 

1000 PIVOTI«PIVCT 

1005 IF(ABS(0ETERM)-R1H030, 1010,1010 
1013 DETERM=DETEPM/r1 


MINV 

10 

^MINV 

20 

MINV 

30 

MINV 

40 

MINV 

50 

MINV 

60 

MINV 

70 

MINV 

80 

MINV 

90 

MINV 

100 

MINV 

110 

MINV 

120 

MINV 

130 

MINV 

140 

MINV 

150 

MINV 

160 

MINV 

170 

MINV 

180 

MINV 

190 

MINV 

200 

MINV 

210 

MINV 

220 

MINV 

230 

MINV 

240 

MINV 

250 

MINV 

260 

MINV 

270 

MINV 

280 

MINV 

290 

MINV 

300 

MINV 

310 

MINV 

320 

MINV 

330 

MINV 

340 

MINV 

350 

MINV 

360 

MINV 

370 

MINV 

380 

MINV 

390 

MINV 

400 

MINV 

410 

MINV 

420 

MINV 

430 

MINV 

440 

MINV 

450 

MINV 

460 

MINV 

470 

MINV 

480 

MINV 

490 

MINV 

500 

MINV 

510 

MINV 

520 

MINV 

530 

MINV 

540 

MINV 

550 

MINV 

560 

MINV 

570 

MINV 

580 

MINV 

590 

MINV 

600 
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1020 


1030 

1040 


1050 

1060 

1070 


1080 


1090 

2000 


2010 

320 


330 

340 

350 

355 

360 

370 


380 

390 

400 

420 

430 

450 

455 

460 

500 

550 


600 

610 

620 

630 

640 

650 

660 

670 

700 

705 

710 

740 


ISCALE*ISCALE+1 

IFCABS( DETERf<l-Rni060f 1020,1020 
0ETERM*DETERM/RI 
ISCALE*ISCAlE4l 
GO TO 1060 

IF( ABS( DETERM )-R2) 1040, 1040,1060 

DETERM«0ETERM*R1 

ISCALE=ISCALE-1 

IF(ABS(0ETEPM J-P2)1050, 1050,1060 

DETcRM*OETERM*Rl 

ISCALE*ISCALE-1 

IF( ABS( PIVOT I J-1^m090, 1070,1070 

PIV0TI*PIV0TI/R1 

ISCALE«ISCAIE+1 

IF(ABS( PIVOTn-RU320,1080, 1080 
PIV0TI*PIV0TI/R1 
ISCALE»ISCALE+1 
GO TO 320 

IF(ABS(PIVOTI )-R2) 2000,2000,320 

PIV0TI»PIV0TI*R1 

ISCALE«ISCALE-1 

IF(ABS(PIVOTI )-R2) 2010, 2010, 320 
PIVOTI*PIVOTI*R1 
ISCALE=1SCALE-1 
DETERM= DETERM* PIVOT I 

DIVIDE PIVOT ROW BY PIVOT ELEMENT 

A( ICOLUM, ICCLUMI^l.O 
DO 350 L*1,N 

A( ICOLUM,LJ = A( ICOLUM, U /PIVOT 
IFCH) 380, 380, 360 
DO 370 L=1,M 

B( ICOLUM,LJ=B( ICOLUM, U /PIVOT 
REDUCE NON-PIVOT ROWS 
DO 550 Ll=l,N 

IF(Ll-ICOLUM) 400, 550, 400 
T=A(L1, ICOLUM) 

A(L1 ,ICOLUM|s.0*0 
DO 450 L=1,N 

A(L1,L)=A(L1,L)-A( ICOLUM,L)*T 
IF<MI 550, 550, 460 
DO 500 L=1,M 

3(L1 ,L)»B(L1,L)-B( ICCLUM,L)*T 
CONTINUE 

INTERCHANGE COLUMNS 

DO 710 1*1, N 
L*N*1-I 

IF (INDEX(L,1I-INDEX(L,2M 630, 710, 
JR0W=INDEX(L,1 ) 

JC0LUM*IN0EX(L,2) 

DO 705 K*1,N 
SWAP=A(K,JROW) 

A(K,JROW)*A(K,JCOLUM) 

A(K,JCOLUM)=SWAP 

CONTINUE 

CONTINUE 

RETURN 

END 


MINV 610 
MINV 620 
MINV 630 
MINV 640 
MINV 650 
MINV 660 
MINV 670 
MINV 680 
MINV 690 
MINV 700 
MINV 710 
MINV 720 
MINV 730 
MINV 740 
MINV 750 
MINV 760 
MINV 770 
MINV 780 
MINV 790 
MINV 800 
MINV 810 
MINV 820 
MINV 830 
MINV 840 
MINV 850 
MINV 860 
MINV 870 
MINV 880 
MINV 890 
MINV 900 
MINV 910 
MINV 920 
MINV 930 
MINV 940 
MINV 950 
MINV 960 
MINV 970 
MINV 980 
MINV 990 
MINVIOOO 
MINVIOIO 
MINV1020 
MINV1030 
MINV1040 
MINV1050 
MINV1060 
MINV1070 
MINV1080 
MINV1090 
MINVllOO 
MINVlllO 
MINV1120 

630 MINV1130 

MINV1140 

MINV1150 

MINV1160 

MINV1170 

MINV1180 

MINV1190 

MINV1200 

MINV1210 

MINV1220 

MINV1230 
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SUBRCUTINE FTLUP ( X » Y, M, N, VARI , V ARO < 

C********OOCUMENT DATE 09-12-69 SUBROUTINE REVISED 07-07-69 


TLUP 


c* 


c* 

C* 


C* 

C* 

c* 


HOOIFICATION OP LIBRARY INTERPOLATION SUBROUTINE 
DIMENSION VAR I (l)fVARO( l),V<3),YY<2> 

DIMENSION IltABI 


FTLUP 


INITIALIZE ALL INTERVAL POINTERS TO -1.0 
DATA <I I( ,J=l»A3)/43*-l/ 

MA*I ABS(M) 


FOR MONOTONICITY 


ASSIGN INTERVAL POINTER FOR GIVEN VARI TABLE 
THE SAME POINTER WILL BE USED ON A GIVEN VARI TABLE EVERY TIME 
LI=MGO(LCCF(VARn n Jf43>n 


i-inu ) 

IF ( I.GE.O) GO TO 10 
IF (N.LT.2) GO TO 10 
C* 

C+MONQTONICITY CHECK 

IF (VARK2J-VARI<1 I ) lfl,3 
C* ERROR IN MCNOTONICITY 

2 K = LOCF (VARKl J) 

PRINT 102f>r,K, CVARIU) , J*=lfN», (VARO< Jl » J*lfN) 

102 FORMAT (IHl,* TABLE BELOW OUT OF ORDER FCR FTLUP AT POSITION 
lrI5f/* X TABLE IS STORED IN LOCATION ♦ » 06 r // f 8G1 5. 0) > 

STOP 

C* MONOTONIC DECREASING 
1 DO 5 J=2fN 

IF <VARI( J)-VARI(J-1))5,2t2 

5 CONTINUE 
GO TO 10 

C* MONOTCNIC INCREASING 

3 DO 6 J=2tN 

IF (VARI(J)-VARI(J-lM2,2t6 

6 CONTINUE 
C* 

C+INTERPOLATICN 

10 IF ( I.LE.O) 1 = 1 
IF (I.GE.N) !=N-1 
IF (N.LE.l) GO TO 8 
IF (MA.NE.O) GO TO 99 
C* ZERO ORDER 
8 Y=VAPD(1I 
GO TO 800 

C* LOCATE I INTERVAL iX< I) .LE . X.LT .X (I ) I 
99 IF nVARK n-X|4(VARI( I+1)-Xn 61,61,40 
C4 IN GIVES DIRECTION FOR SEARCH OF INTERVALS 

40 IN=SIGN(1.0,(VARI( I ♦! I -VARI ( 1 1 )♦ (X-VARI ( I )) ) 

C* IF X OUTSIDE ENDPOINTS, EXTRAPOLATE FROM END INTERVAL 

41 IF ni + IN).LE.O) GO TO 61 
IF (d+IN) .GE.NJ GO TO 61 
1=I+IN 

IF UVARK n-X)*(VARin+l)-Xn 61,61,41 
IF (MA.EQ.2) GO TO 200 


61 
C* 
C^FIRST 


ORDER 


TLUP 
TLUP 
TLUP 
TLUP 
CHECKTLUP 
TLUP 
TLUP 
TLUP 
TLUP 
TLUP 
TLUP 
TLUP 
TLUP 
TLUP 
TLUP 
TLUP 
TLUP 
TLUP 
TLUP 
TLUP 
TLUP 
TLUP 
TLUP 
TLUP 
TLUP 
TLUP 
TLUP 
TLUP 
TLUP 
TLUP 
TLUP 
TLUP 
TLUP 
TLUP 
TLUP 
TLUP 
TLUP 
TLUP 
TLUP 
TLUP 
TLUP 
TLUP 
TLUP 
TLUP 
TLUP 
TLUP 
TLUP 
TLUP 
TLUP 
TLUP 
TLUP 
TLUP 
TLUP 


Y=(VARO( n*(VAPl(I*n-X)-VAPO( l4-l)*(VARI(n-Xn/<VARin+l)-VART( DTLUP 


1 ) 

GO TO 800 
C* 

C+ScCONO ORDER 


TLUP 

TLUP 

TLUP 

TLUP 


10 

20 

30 

40 

50 

60 

70 

80 

90 

100 

110 

120 

130 

140 

150 

160 

170 

180 

190 

200 

210 

220 

230 

240 

250 

260 

270 

280 

290 

300 

310 

320 

330 

340 

350 

360 

370 

380 

390 

400 

410 

420 

430 

440 

450 

460 

470 

480 

490 

500 

510 

520 

530 

540 

550 

560 

570 

580 

590 

600 
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200 IF (N.EQ.2I GO TO 2 

IF ( I.EQ.(N-l)) GO TO 209 
IF ( I. EG. II GO TO 201 
C* PICK THIRD POINT 

SK= VARI ( I*1I-VARI ( 1 1 

IF ( (SKfl K-VAPI ( 1-1 ) 1 1 .LT. ( SK* ( VARI ( H-2)-Xl II 


GO TO 209 


201 

209 

702 


800 


L = I 

GO TO 702 
L = I-1 

V(11=VAR1(L1-X 

VI2I=VARHL*1I-X 

YY(1I = ( VAR0(LI*V(2 l-VARD(L*l I*V( II I / ( VAR I ( L+1 l-VARH L I I 

YY(2I={ VARDIL*11*VI 31-VARO(L*2I*V(2n/(VARI(L + 2l-VARI (L+lII 

Y= (YY( I l♦V(3I-YY^2 I*V( 1 ) I/( VARI ( L+2 l-VARI< LI I 


II(LII = I. 

RETURN 

ENO 


TLUP 610 
TLUP 620 
TLUP 630 
TLUP 640 
TLUP 650 
TLUP 660 
TLUP 670 
TLUP 680 
TLUP 690 
TLUP 700 
TLUP 710 
TLUP 720 
TLUP 730 
TLUP 740 
TLUP 750 
TLUP 760 
TLUP 770 
TLUP 780 
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Figure 3.- This detailed sketch of a chordwise row of horseshoe vortices illustrates the 
velocities and circulations used to compute lift and pitching moment on the elemental 
panels of a wing with dihedral. Note that the velocity terms and circulations which 
are shown with each horseshoe vortex are different. (See Part m, Section 1 for 
discussion.) 
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Figure 4.- Spanwise bound vortex filament at an arbitrary orientation in the flow. 
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Lift computed on trailing vortex filament 
Data from output listing 

Lift from span wise vortex filament 
Lift from trailing vortex filament 
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Figure 7.- Concluded. 









Figure 8.- Concluded. 








Figure 9.- Concluded. 
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